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ABSTRACT 


A  solution  has  been  obtained,  to  the  magnetogasdynamie 
equations  of  motion  including  the  effects  of  compressibility, 
viscosity,  thermal  conductivity  and  crossed  electric  and 
magnetic  fields.  The  solution  obtained  is  valid  for  an 
acceleration  flow  in  which  a  separation  of  variable  form 
of  solution  is  utilized.  The  constraints  under  which  a 
separable  solution  is  valid  are  clearly  delineated.  The 
solution  is  an  exact  solution  in  the  sense  that  no  terms 
have  been  dropped  from  the  governing  equations  and  boundary 
layer  approximations  have  not  been  employed. 

The  results  of  the  investigation  were  presented  in  the 
form  of  velocity  and  temperature  profiles  for  several  values 
of  the  electromagnetic  parameters  Ha  and  $.  Although  the 
mathematical  model  used  in  the  present  investigation  is  too 
crude  to  be  utilized  for  detailed  design  of  a  practical  MHD 
accelerator,  it  is  believed,  nevertheless,  that  the  model 
displays  certain  qualitative  features  of  the  flow  which  are 
likely  to  be  encountered  in  a  practical  device. 


iii 


AEDC-T  R-67-80 


TABLE  OF  CONTENTS 

Page 

ABSTRACT .  iii 

NOMENCLATURE  .  xi 

I.  INTRODUCTION . . .  1 

II.  PREVIOUS  WORK 

2.1  One  Dimensional  Studies  .  4 

2.2  Incompressible  Studies  . 5 

2.3  Viscous  Magnetogasdynamic  Studies  .  6 

III.  FORMULATION  OF  THE  PROBLEM  .  9 

IV.  TECHNIQUE  OF  SOLUTION 

4.1  Solution  Difficulties  of  the  Direct 

Boundary  Value  Problem  . 23 

4.2  The  Initial  Value  Problem  . .  28 

V.  DISCUSSION  OF  RESULTS  .  37 

VI.  VARIABLE  PROPERTIES 

6.1  Formulation  of  the  Equations  .  44 

6.2  Discussion  of  Results  .  47 

VII.  SUMMARY  .  49 

VIII.  REFERENCES  .  51 


v 


AEDC-TR-67-80 


LIST  OF  ILLUSTRATIONS 

Figure 

1.  MHD  Powered  Wind  Tunnel 

2.  Channel  Configuration 

3.  Relationship  Among  Parameters  a,  Ha  and  <£> 

4.  Relationship  Among  Parameters  g(a),  Ha,  O, 
and  a 

5.  Region  of  Acceptable  Solutions  for  Ha  =  200, 

0  =  1.6 

6.  Velocity  and  Temperature  Profiles  for  Ha  =  100, 

<£  =  1.2 

7.  Velocity  and  Temperature  Profiles  for  Ha  =  100, 

<t>  =  1.6 

8.  Velocity  and  Temperature  Profiles  for  Ha  =  100, 

3>  =  2.0 

9.  Velocity  and  Temperature  Profiles  for  Ha  -  150, 

$  -  1.2 

10.  Velocity  and  Temperature  Profiles  for  Ha  =  150, 

$  =  1.6 

11.  Velocity  and  Temperature  Profiles  for  Ha  =  150, 

$  =  2.0 

12.  Velocity  and  Temperature  Profiles  for  Ha  =  200, 

$  =  1.2 

13.  Velocity  and  Temperature  Profiles  for  Ha  =  200, 

0  =r  1.6 

14.  Velocity  and  Temperature  Profiles  for  Ha  =  200, 

$  =  2,0 

15.  Variation  of  the  Velocity  Gradient  at  the  Wall 
for  a  =  1.5°,  g(a)  =  0.2 

16.  Variation  of  the  Temperature  Gradient  at  the 
Wall  for  a  =  1.5°,  g(a)  =  0.2 

17.  Velocity  and  Temperature  Profiles  for  Ha  “  100, 
<J>  =  1.6 

18.  Velocity  and  Temperature  Profiles  for  Ha  =  150, 

ct  =  1.2 


vi 


19. 

Velocity  and 
$  =  1.6 

20. 

Velocity  and 
^  -  1.6,  o  ~ 

21. 

Velocity  and 
<2>  =  1.2,  o  ^ 

22. 

Velocity  and 
0=1.6,  a 

23. 

Velocity  and 

0  =  1.6,  0  'v. 

24. 

Velocity  and 
0  =  1,2,  0  '■>-< 

25. 

Velocity  and 
0  =  1.6,  0  ^ 

26. 

Variation  of 
0  =  1.2 

27. 

Variation  of 
0  =  1.2 

28. 

Variation  of 
0  =  1.6 

29. 

Variation  of 
0  =  1.6 

30. 

Variation  of 
0  =  2.0 

31. 

Variation  of 
0  =  2.0 

32. 

Variation  of 
0=1.2 

33. 

Variation  of 
0  =  1.2 

34. 

Variation  of 
0  =  1.6 

35. 

Variation  of 
0  =  1.6 

36. 

Variation  of 
0  =  2.0 

A  E  DC-T  R-67-80 

Temperature 

Profiles 

for 

Ha  = 

150, 

Tenjgerature 

Profiles 

for 

Ha  = 

100, 

Temperature 

T3/2 

Profiles 

for 

Ha  = 

150, 

Temperature 

T3/2 

Profiles 

for 

Ha  = 

150, 

Temperature 
T3/2,  p  ^  K 

Profiles 
^  T1/2 

for 

Ha  - 

100, 

Temperature 
T3/  ,  [l  ^  K 

Profiles 

^  T1/2 

for 

Ha  = 

150, 

Temperature 

T3/2,  jj,  ^  K 

Profiles 
^  T1?2 

for 

Ha  = 

150, 

a  with  Re  and  M  for  Ha  =  100, 


g(a)  with  Re  and  M  for  Ha  -  100, 
a  with  Re  and  M  for  Ha  =  100, 
g(a)  with  Re  and  M  for  Ha  =  100, 
a  with  Re  and  M  for  Ha  =  100, 
g(a)  with  Re  and  M  for  Ha  =  100, 
a  with  Re  and  M  for  Ha  =  150, 
g(a)  with  Re  and  M  for  Ha  =  150, 
a  with  Re  and  M  for  Ha  =  150, 
g(a)  with  Re  and  M  for  Ha  =  150, 
a  with  Re  and  M  for  Ha  =  150, 


Vll 


AEDC-TR-67-B0 


37. 

Variation 

of 

g<a) 

with 

4> 

=  2.0 

38. 

Variation 

of 

a  with  Re 

4> 

=  1.2 

39. 

Variation 

of 

g(a) 

with 

4> 

=  1.2 

40. 

Variation 

of 

a  with  Re 

4> 

=  1.6 

41. 

Variation 

of 

g(a) 

with 

4> 

=  1.6 

42. 

Variation 

of 

a  with  Re 

4> 

=  2.0 

43. 

Variation 

of 

g(a) 

with 

4> 

=  2.0 

44. 

Variation 

of 

a  with  Re 

4> 

-  1.6,  , 

a  '"v 

T3/2 

45. 

Variation 

of 

g(a) 

with 

4> 

=  1.6, 

0  nj 

T3/2 

46. 

Variation 

Of 

a  with  Re 

4> 

=  1.2, 

0  <tl 

T3/2 

47. 

Variation 

Of 

g(a) 

with 

4> 

=  1.2, 

0  'v. 

T3/2 

48. 

Variation 

of 

a  with  Re 

$ 

=  1.6, 

cr  'vi 

T3/2 

49. 

Variation 

of 

g  (a) 

with 

4> 

=  1.6, 

O'  'v 

T3/2 

50. 

Variation 

of 

a  with  Re 

4> 

=  1.6, 

0  'v 

T  3/?, 

n  ~ 

51. 

Variation 

of 

g(a) 

with 

4> 

=  1.6, 

0  'V. 

T3/2, 

Ul 

to 

• 

Variation 

of 

a  with  Re 

0 

=  1.2, 

0  'Vi 

T3/2 , 

53. 

Variation 

of 

g(a) 

with 

4> 

-  1.2, 

a  <\i 

T3/2, 

|X  ~ 

54. 

Variation 

of 

a  with  Re 

4> 

-  1.6, 

0 

T3/2, 

|i  ^ 

viii 


Re  and  M  for  Ha  =  150, 

and  M  for  Ha  =  200, 

Re  and  M  for  Ha  =  200, 

and  M  for  Ha  =  200, 

Re  and  M  for  Ha  =  200, 

and  M  for  Ha  =  200, 

Re  and  M  for  Ha  -  200, 

and  M  for  Ha  =  100, 

Re,  M  for  Ha  =  100, 

and  M  for  Ha  =  150, 

Re  and  M  for  Ha  =  150, 

and  M  for  Ha  =  150, 

Re  and  M  for  Ha  =  150, 

and  M  for  Ha  =  100, 
k  -  T1/2 

Re  and  M  for  Ha  =  100, 

k  ~  T1/2 

and  M  for  Ha  =  150, 

K  ~  T1/2 

Re  and  M  for  Ha  =  150, 

K  -v,  T1/2 

and  M  for  Ha  =  150, 

K  *  T1/2 


AEDC-TR-67-80 


55. 


Variation  of  g(a)  with  Re  and  M  for  Ha 
$  =  1.6,  o  'v  T3/2 ,  |i  -v  (c  -v  T1/2 


150, 


AEDC-TR-67-8Q 


I 

II 

III 


LIST  OF  TABLES 

Limiting  Values  of  a,  M  and  Re 

Summary  of  Results  for  Constant  Properties 

Comparison  of  f^a)  and  gf(a)  for  Constant  and 
Variable  Properties 


x 


AEDC-TR-67-B0 


B 

cp 

d 

-> 

E 

f(0) 

g(0) 

h 

Ha 

Ha* 

-> 

J 

M 

P 

Pr 


r 

R 

Re 

Re* 

RM 

T 

T(0) 

-> 

U 

u  (0) 

z 


a 

7 


NOMENCLATURE 

Magnetic  induction  (B^,  B^,  Bz) 
Specific  heat  capacity 
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Electric  field  vector  (E  ,  E„,  E  ) 

X  w  Z 

Dimensionless  velocity,  Eq.  (42) 
Dimensionless  temperature,  Eq.  (42) 
Enthalpy 

Hartmann  number,  Eq.  (44) 

Hartmann  number  based  on  dQ 

Current  density  vector  (J  ,  Ja,  J  ) 

r  t7  z 

Mach  number,  Eq.  (44) 

Static  pressure 
Prandtl  number,  Eq .  (44) 

Coordinate  direction 

Gas  constant 

Reynolds  number,  Eq.  (44) 

Reynolds  number  based  on  dQ 
Dimensionless  parameter,  Eq,  (60) 
Temperature 

Dimensional  temperature  function 
Velocity  vector  (U  ,  Ua,  U) 
Dimensional  velocity  function 
Coordinate  direction 
Divergence  angle 
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ri  Normalized  coordinate,  0/a 

0  Coordinate  direction 

k  Thermal  conductivity 

X  Mean  free  path 

p,  Viscosity 

(ie  Magnetic  permeability 

it  Dimensionless  parameter,  Eq.  (60) 

p  Density 

a  Electrical  conductivity 

4>  Loading  parameter,  Eq.  (44) 

a)  Exponent,  Eq.  (110) 


SUBSCRIPTS 

c  Characteristic  or  reference  values 

o  Values  at  the  channel  entrance,  r  =  rQ 

w  Values  at  the  wall  9  =  a 
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I 

INTRODUCTION 


For  several  years  there  has  been  considerable  interest 
in  the  electromagnetic  acceleration  of  ionized  gases  for 
space  propulsion  purposes  and  for  application  to  high  speed 
wind  tunnels.  Many  types  of  devices  have  been  proposed  to 
accomplish  this  acceleration  and  a  number  of  prototypes 
have  been  built  that  demonstrate  technical  feasibility. 

One  of  the  more  promising  of  these  devices  is  the  steady 
flow,  crossed  field  or  Faraday  accelerator.  There  is 
particular  interest  in  this  type  of  accelerator  for  wind 
tunnel  applications. 

One  factor  that  limits  conventional  hypervelocity 
flight  simulation  techniques  is  the  necessity  of  confining 
a  high  temperature  gas  at  high  pressure  prior  to  expanding 
in  a  supersonic  nozzle.  By  directly  adding  kinetic  energy 
to  the  flow,  a  magnetogasdynamic  accelerator  would,  to  a 
large  extent,  alleviate  this  containment  problem.  Ring 
(Ref.  1)  has  stated,  "In  a  sense,  MHD  techniques  offer  the 
first  promise  of  a  quantum  increase  in  aerodynamic  test 
capabilities  since  the  introduction  of  shock- tube,  shock- 
tunnel  techniques  a  number  of  years  ago." 

A  schematic  drawing  of  a  crossed  field  accelerator  to 
be  used  for  the  simulation  of  high  speed  flight  is  shown 
in  Fig.  1.  In  such  an  arrangement  air  would  be  preheated, 
seeded,  and  ionized  in  an  arc  heater  and  then  passed  through 
a  low  Mach  number  supersonic  nozzle  into  the  accelerating 
channel.  In  the  acceleratoi’  itself,  energy  would  be  added 
to  the  ionized  gas  or  plasma  by  the  crossed  electric  and 
magnetic  fields.  A  large  portion  of  this  added  energy 
would  be  In  the  form  of  directed  kinetic  energy,  the  re¬ 
mainder  being  random  thermal  energy  increasing  the  enthalpy 
of  the  plasma.  A  post-accelerator  nozzle  would  expand  the 
gas  to  the  desired  static  conditions  and  further  increase 
the  Mach  number. 

It  is  clear  that  the  accelerator  itself  is  the  heart 
of  such  a  flight  simulation  system  and  that  it  is  important 
that  the  details  of  the  flow  in  this  device  be  well  under¬ 
stood  . 

A  considerable  amount  of  theoretical  work  has  been 
done  on  the  crossed  field  accelerator  and  in  the  next  section 
a  brief  survey  of  the  more  pertinent  of  these  analyses  is 
presented.  The  majority  of  these  investigations  have  treated 


1 


AE  DC-T  R  '67»80 


an  inviscid,  non-heat  conducting  gas,  the  tacit  assumption 
being  made  that  the  viscous  and  heat  transfer  effects  are 
confined  to  a  thin  boundary  layer  adjacent  to  the  channel 
wall.  However,  in  wind  tunnel  accelerators  where  relatively 
large  length  to  diameter  ratios  will  be  necessary  to  produce 
the  required  acceleration,  it  is  likely  that  a  significant 
portion  of  the  flow  will  be  fully  developed  in  the  sense 
that  the  boundary  layers  have  merged.  Thus,  the  influences 
of  viscosity  and  thermal  conductivity  would  extend  throughout 
the  flow  field.  This  would  especially  be  true  for  accelerators 
operating  at  low  to  moderate  Reynolds  number  i.e.  low  density. 
It  is  well  known  (Ref.  2)  that  under  these  conditions  the 
flow  must  be  treated  as  fully  viscous  and  that  boundary 
layer  theory  is  no  longer  applicable.  Under  these  circum¬ 
stances,  the  results  of  one-dimensional  analyses  cannot  be 
expected  to  accurately  describe  the  actual  flow. 

It  would  appear  then  that  there  is  a  need  for  a  study 
of  the  crossed  field  accelerator  including  the  effects  of 
compressibility,  viscosity,  and  thermal  conductivity.  A 
general  study  of  such  a  problem  would  be  a  formidable  under¬ 
taking  indeed.  In  this  investigation  attention  is  focused 
on  a  special  type  of  accelerator  flow  for  which  it  is 
possible  to  obtain  a  solution  to  the  governing  equations 
including  the  aforementioned  effects. 

The  flow  considered  is  a  two-dimensional  gas  flow 
between  diverging  plane  walls,  the  pressure  and  Mach  number 
being  constant  down  the  channel.  For  this  model  it  is 
possible  to  separate  the  axial  and  cross  channel  dependence 
of  velocity  and  temperature.  Performing  this  separation 
results  in  a  pair  of  coupled,  non-linear,  ordinary  differential 
equations  for  the  cross  channel  variations.  This  system  of 
equations  is  extremely  complex  and  must  be  integrated  by 
numerical  methods.  Moreover  it  is  shown  that  the  problem 
is  generally  over  specified  and  that  solutions  cannot  be 
obtained  for  arbitrary  values  of  the  parameters.  This 
difficulty  is  circumvented  by  converting  the  boundary  value 
problem  to  an  initial  value  problem  and  allowing  the  wall 
divergence  angle  and  wall  temperature  to  be  determined  by 
the  integration. 

In  order  to  accomplish  the  separation  of  variables 
mentioned  above  it  is  necessary  to  assume  that  the  viscosity, 
thermal  conductivity,  and  electrical  conductivity  do  not 
change  down  the  channel.  Consequently,  in  the  main  portion 
of  this  report,  it  is  assumed  that  these  properties  are 
constant  throughout  the  flow  field.  This  assumption  rep¬ 
resents  a  rather  severe  restriction  bfecause  the  viscosity 
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and  electrical  conductivity  are  strong  functions  of  tem¬ 
perature.  The  latter  section  of  this  report  is  concerned 
with  an  attempt  to  assess  the  effect  of  these  property 
variations  across  the  channel  on  the  velocity  and  temperature 
profiles . 
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II 

PREVIOUS  WORK 


2.1  ONE-DIMENSIONAL  STUDIES 

The  one-dimensional  magnetogasdynamic  equations  for  an 
inviscid  non-heat  conducting  gas  were  presented  by  Resler 
and  Sears  (Ref,  3),  and  many  of  the  interesting  properties 
of  these  equations  were  pointed  out.  In  particular,  the 
possibility  of  accelerating  a  gas,  flowing  in  a  channel,  to 
very  high  speeds  was  noted.  In  another  paper  (Ref.  4)  the 
same  authors  discussed  in  some  detail  constant  area  magneto- 
gasdynamic  channel  flow  listing  eight  possible  flow  situations 
that  might  occur,  depending  on  the  initial  velocity  and  initial 
Mach  number. 

In  the  succeeding  years,  a  number  of  authors  have  analyzed 
quasi -one-dimensional  channel  flow  problems  with  crossed 
electric  and  magnetic  fields.  Matthews  (Ref.  5),  for  example, 
has  presented  a  number  of  solutions  for  the  steady  flow  crossed 
field  accelerator.  Following  the  approach  of  Resler  and  Sears, 
Matthews  assumes  the  electrical  conductivity  to  be  a  constant 
and  ignores  the  effect  of  the  induced  magnetic  field.  The 
latter  assumption  permits  the  gas  dynamical  equations  to  be 
treated  independently  of  Maxwell's  electromagnetic  equations. 

A  further  consequence  of  this  assumption  is  that  the  variation 
of  two  of  the  dependent  variables  must  be  specified  before  a 
solution  can  be  obtained. 

Drake  (Ref.  6)  studied  the  isothermal  acceleration  of 
a  plasma  with  the  aim  of  optimizing  the  channel  design.  Using 
a  variational  technique  he  determined  the  channel  shape  that 
would  yield  a  minimum  accelerator  length  for  a  given  velocity 
change. 

Krupka  and  Kezios  (Ref.  7)  analyzed  a  constant  area 
accelerator  for  constant  applied  electric  and  magnetic  fields. 
In  their  study  the  induced  field  in  the  direction  of  the 
applied  field  was  not  neglected  and  the  electrical  conductivity 
was  taken  to  be  a  function  of  temperature.  The  results  of 
Krupka  and  Kezios  indicate  that  the  Mach  number  increase 
during  acceleration  is  limited  and  that  the  flow  cannot  be 
accelerated  through  the  sonic  velocity. 

Two  investigations  of  inviscid  magnetogasdynamic  flow 
are  of  particular  interest  in  connection  with  the  present 
work.  These  are  the  papers  by  Hains  (Ref.  8)  and  Podolsky 
and  Borman  (Ref.  9).  Hains  examined  several  cases  of 
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magnetogasdynamic  source  flow  and  obtained  a  particularly 
simple  solution  for  the  case  of  crossed  electric  and  mag¬ 
netic  fields.  The  electric  field  was  assumed  constant  and 
the  magnetic  field  was  assumed  to  vary  inversely  with 
radial  distance  from  the  source.  The  electrical  con¬ 
ductivity  was  assumed  constant  and  induced  field  effects 
and  Hall  effects  neglected.  Hains  found  that  the  flow 
proceeds  outward  from  the  source  at  constant  pressure  and 
constant  Mach  number,  the  velocity  increasing  linearly 
with  radial  distance  and  the  temperature  increasing  as  the 
square  of  the  radial  distance. 

In  a  related  investigation,  Podolsky  and  Borman  (Ref. 

9)  considered  precisely  the  geometry  that  is  dealt  with 
in  the  present  study.  The  governing  equations  for  a  crossed 
field  accelerator  were  written  in  cylindrical  coordinates, 
neglecting  changes  in  the  z  and  9  directions.  These  equations 
were  then  numerically  integrated  for  specified  entrance 
conditions  at  r  .  For  large  values  of  r  they  obtained  a 
solution  of  thexorm 


u  ^  r 

p  =  constant 
T  ~  r2 


From  the  solution  graph  presented  it  appears  that  this 
pattern  is  established  fairly  quickly  and  is  valid  for 
values  of  r/rp  >3.  In  the  entrance  region  where  1  <  r/rQ  <  3 
the  velocity  is  still  almost  linear  although  the  pressure 
drops  rapidly.  Podolsky  and  Borman  also  treat  this  problem 
for  a  lithium  vapor  propellant  using  values  for  the  con¬ 
ductivity  calculated  from  the  degree  of  ionization  at 
equilibrium. 

Sherman  (Ref.  10)  has  extended  the  work  of  Podolsky 
and  Borman  to  include  Hall  currents  and  propellants  that 
are  mixtures  of  inert  gas  and  ionized  seed  gas. 

2.2  INCOMPRESSIBLE  STUDIES 

A  vast  number  of  articles  have  been  published  that 
deal  with  magnetohydrodynamic  channel  flows,  and  only  a  few 
of  the  more  pertinent  will  be  discussed  here.  No  survey  of 
the  literature  in  this  area  would  be  complete,  however, 
without  mentioning  the  original  paper  in  the  field  by  Hartmann 
(Ref.  11),  which  deals  with  fully  developed  laminar  flow 
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between  parallel,  non-conducting  walls.  A  discussion  of 
Hartmann's  work  and  several  other  interesting  magnetohydro¬ 
dynamic  channel  flows  is  contained  in  a  recent  book  by 
Sutton  and  Sherman  (Ref.  12). 

In  an  incompressible  study  somewhat  similar  to  the 
present  work,  Axford  (Ref.  13)  examined  the  magnetohydro¬ 
dynamic  flow  between  non-parallel  plane  walls,  an  extension 
of  the  classical  Jef f rey-Hamel  problem.  The  channel  shape 
and  magnetic  field  distribution  were  the  same  as  those 
considered  in  this  investigation,  but  the  applied  electric 
field  was  taken  to  be  zero,  corresponding  to  a  short  cir¬ 
cuited  MHD  generator.  Writing  the  governing  equations  in 
cylindrical  coordinates,  Axford  was  able  to  separate  the 
axial  and  cross  channel  dependence.  By  assuming  small 
magnetic  Reynolds  number,  i.e.  neglecting  the  induced  field, 
he  was  able  to  obtain  a  closed  form  solution  for  the  cross 
channel  variation.  The  emphasis  in  this  study  was  on  the 
mathematical  analysis  and  only  a  limited  number  of  numerical 
results  were  presented.  It  was  reported  that  for  divergent 
flow  the  tendency  toward  separation  and  back  flow  at  high 
Reynolds  number  is  opposed  by  the  transverse  magnetic  field. 
Velocity  profiles  presented  for  convergent  flow  show  the 
characteristic  flattening  for  increasing  magnetic  field 
strength. 

Vatazhin  (Ref.  14)  examined  some  limiting  cases  of 
the  divergent  channel  problem  treated  by  Axford  and  obtained 
solutions  in  terms  of  elementary  functions,  Vatazhin 's 
results  support  the  conclusions  made  by  Axford. 

Heywood  (Ref.  15)  investigated  the  magnetohydrodynamic 
flow  between  parallel  plates  assuming  an  electrical  con¬ 
ductivity-temperature  relation  of  the  form  c  ~  T^.  He 
found  very  little  effect  on  the  temperature  profile  due 
to  the  conductivity  variation  but  a  noticeable  alteration 
of  the  velocity  profile.  He  noted,  however  that  for  gases 
flowing  at  high  speed  the  effect  on  the  temperature  dis¬ 
tribution  could  be  significant. 

2.3  VISCOUS  MAGMETOGASDYNAMIC  STUDIES 

In  an  early  paper,  Bleviss  (Ref.  16)  studied  magneto- 
gasdynamic  Couette  flow  in  an  effort  to  deduce  information 
about  the  magnetogasdynamic  boundary  layer.  Hall  effects 
were  not  included  in  this  study  but  the  variation  of  trans¬ 
port  properties  was  taken  into  account.  Bleviss'  results 
show  a  rather  severe  distortion  of  the  velocity  and  temper¬ 
ature  distributions  with  increasing  magnetic  field  strength. 
This  study  was  directed  toward  aerodynamic  boundary  layers 
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on  space  vehicles  reentering  the  atmosphere  rather  than 
boundary  layers  in  channels.  Consequently  the  effect  of 
an  applied  electric  field  was  not  considered. 

Martin  (Ref.  17)  studied  the  flow  of  a  compressible, 
viscous  fluid  in  a  pipe  of  circular  cross  section;  the 
fluid  being  accelerated  by  an  axial  body  force.  His  work 
dealt  with  both  the  entrance  flow  and  the  fully  developed 
flow  which  was  defined  as  a  non-accelerating  parallel  flow. 

A  body  force  is  required  in  order  for  the  fully  developed 
flow  to  exist.  By  assuming  that  the  viscosity,  thermal 
and  electrical  conductivities  were  constant  and  that  the 
flow  was  unidirectional  throughout,  Martin  was  able  to 
solve  the  linearized  momentum  equation  independently  of 
the  energy  equation. 

Although  the  applicability  of  this  study  to  magneto- 
gasdynamic  channel  flows  was  stressed,  it  is  doubtful  that 
the  body  forces  assumed  could  actually  be  produced  by  the 
interaction  of  electric  and  magnetic  fields  in  circular 
pipes . 

Some  very  interesting  results  were  obtained  by  Hale 
and  Kerrebrock  (Ref.  18)  in  a  study  of  insulator  boundary 
layers  in  crossed  field  accelerators.  Two  models  for  the 
electrical  conductivity  variation  were  employed;  an  equi¬ 
librium  model  where  the  degree  of  ionization  was  determined 
by  the  local  gas  temperature  and  pressure,  and  a  non¬ 
equilibrium  model  in  which  the  electrical  conductivity  is 
strongly  coupled  to  the  electric  field  strength. 

It  was  found  that  similar  solutions  could  be  obtained 
only  for  very  restricted  conditions  and  local  similarity 
was  assumed. 

The  boundary  layer  profiles  for  the  non-equilibrium 
case  show  some  very  unusual  characteristics  and  were  ex¬ 
tremely  sensitive  to  the  free  stream  Mach  number  and  elec¬ 
tromagnetic  loading  parameter.  For  certain  values  of  the 
parameters,  the  velocity  and  temperature  profiles  show  very 
pronounced  bulges  near  the  wall.  The  variation  of  the 
profiles  for  the  equilibrium  case  is  less  pronounced. 

Sonnerup  (Ref.  19)  investigated  the  fully  developed 
flow  of  a  compressible,  electrically  conducting  fluid  in 
slowly  diverging  channels.  Sonnerup' s  work  differs  from 
the  present  study  in  that  he  was  concerned  with  flow  in  a 
generator  and  the  confining  walls  were  electrode  surfaces 
rather  than  insulator  surfaces.  The  most  interesting  feature 
of  this  work  is  that  a  closed  form  solution  was  obtained. 
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Sonnerup  simplified  the  governing  equations  via  boundary 
layer  type  approximations  in  a  manner  similar  to  Williams 
(Ref.  2).  By  choice  of  the  proper  wall  geometry  the  pressure 
gradient  and  the  electromagnetic  body  force  were  made  to 
balance  each  other.  This,  in  turn,  permitted  the  energy 
equation  and  the  momentum  equation  to  be  uncoupled  and 
analyzed  separately.  The  velocity  profiles  were  found  to 
be  self-similar  but  the  temperature  profiles  were  generally 
non-similar . 

Although  the  emphasis  of  SonnerupJs  work  was  on  the 
mathematical  analysis  some  numerical  results  were  presented. 
The  velocity  and  temperature  profiles  given  exhibit  no 
unusual  characteristics. 
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III 

FORMULATION  OF  THE  PROBLEM 


Consider  the  steady,  two-dimensional,  laminar  flow  of 
an  electrically  conducting  gas  between  diverging  plane  walls 
as  shown  in  Fig.  2,  The  problem  is  best  described  in  cylin¬ 
drical  polar  coordinates. 

The  flow  may  either  be  considered  as  emerging  from  a 
line  source  at  r  =  0  or  as  the  fully  developed  portion  of 
an  actual  channel  flow.  The  rays  9  =  +  a  are  electrically 
insulated  walls  and  the  magnetic  field  is  applied  in  the  9 
direction.  An  electric  field  is  applied  in  the  negative  z 
direction,  i.e.  into  the  paper. 

The  governing  equations  of  motion  for  a  single  fluid 
plasma  are  expressed  in  cylindrical  polar  coordinates  in 
Eqs.  (1)  through  (5). 


Continuity  equation: 


dp  .  1  s<pV>  .  1  a<PV  .  3(P*V  „ 

-  +  —  ■  ■  ■  +  —  *  +  -  2=2  V 


dt  r  or 


r  be 


dz 


(1) 


r  momentum  equation: 


DU  V  r?  +  *  ^  *  dU 

p[—  -  JL]  =  (JxB)  -  ^  +  -2-  [|i (2  — 
Dt  r  dr  dr  dr 


2  v 


U)] 


,  N  -i  dU  dUfl  Ua  ■* 

1  d  r  /I  4.  — ®  °  r  f  r  ■  z 


+  ±  —  [n(- 

r  d@  t  bO  dr 


dU  dU 

)]  +f" 

r  dz  oz  or 


+  at  (5e  _  1  ^ _  % 

r  dr  r  d0  r 


(2) 
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0  momentum  equation: 


P[^S  +  ^2]  =  + 

Dt  r  r  50  r  50  r  5 0 


2U  „  ->  >  5Ufl  ,  5U  N  ,  ^11 

+  —  -  -  V  •  U)]  +  -  +  -  — ^)]  +  -i- [n(-  — 

r  3  5z  5z  r  50  5r  r  50 


dUn  2|1  ,  5U  5Ufl  Uft 

-  —2)3  +  (I  _I  +  -1  -  _£) 


d  r 


r  r  50  5r 


(3) 


z  momentum  equation: 


DUz 

p  — -  =  (JxB) 


z  _  d£ 


Dt 


-v 

+  ^-[n(2  — ^ 
5z  5z  5z 


o  -» 

-  V  .  U) 


0. 


,  i  5u  5u  du„  5U 

+  —  -^-[tirC—  +  — )]  +  -  —  +  — )] 

r  5r  5z  5r  r  50  r  50  5z 


(4) 


Energy  equation: 


P  a  =  Se.  +  2%  I  A  („  21,  +  i  A(k  il,  +  A.u  3t> 

Dt  Dt  c  r  5r  dr  r2  50  50  5z  5z 


+  u 


r  5U  2  -i  5Ufl  U 

|2  ( — — )  +2(i— S  + 

[  5r  r  50  r 


a  SU  2  i  5U  5U0 

•)  +  2  ( — —)2  +  (i  — ^  +  _^)2 


a  z 


r  50  5z 


SU  5U  a  !  SU  due  Ua  ■ 

+  (—  +  — )2  +  (i  —  +  —  +  — ) 


5z  5  r 


r  50  5r 


-  -  |x(V  •  U)2  (5) 
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_D_ 

Dt 


A 

Sr  r 


U. 


_S_ 

Sz 


V 


u 


!  S(rUr) 
r  Sr 


r  S  0 


The  necessary  equations  for  the  electromagnetic  variables 
are : 


Ampere's  laws 


V  x  B  =  |ieJ 


(6) 


Faraday’s  law: 


V  x  E 


SB 

St 


(7) 


Ohm ’ s  law : 


J  =  a(E  +  U  x  B) 


(8) 


The  equations  as  written  incorporate  the  so-called 
magnetohydrodynamic  approximations,  i.e.  the  dielectric 
displacement  current  is  neglected  in  Eq.  (6)  and  the  gas 
is  assumed  to  be  electrically  neutral.  The  latter  assumption 
permits  deletion  of  the  convection  current  term  in  Ohm's  law 
and  a  body  force  term  due  to  the  electric  field.  These 
approximations  are  discussed  by  Sutton  and  Sherman  (Ref.  12) 
and  are  standard  magnetof luiddynamic  assumptions. 

Ohm’s  law  is  written  in  its  simplified  form  neglecting 
Hall  currents  and  ion  slip  effects. 

In  addition  to  the  equations  listed  above  we  have  an 
equation  of  state: 


p  =  p(p,T) 


(9) 
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a  caloric  equation  of  state: 


h  =  h  (p ,  T)  (10) 


and  relations  for  the  transport  properties: 


li  =  «  0  =  cr(p,T)  ,  k  =  <  (p ,  T)  (11) 


Eqs.  (1)  through  (10)  ar§  sixteen  equations  for  the 
sixteen  unknowns  p,  p,  T,  h,  U,  B,  £,  and  J.  With  appropriate 
boundary  conditions  they  represent,  theoretically,  a  deter¬ 
minate  set.  These  equations  will  now  be  reduced  to  a  form 
suitable  for  the  particular  problem  being  investigated. 

As  was  stated  previously,  the  flow  is  assumed  to  be 
steady  and  two  dimensional,  i.e.  d/dt  =  d/dz  =  0.  Also,  the 
transverse  velocity  components  U@  and  Uz  are  assumed  to  be 
zero.  It  is  further  assumed  that  the  induced  magnetic  field 
is  negligibly  small  compared  to  the  applied  magnetic  field. 
This  last  assumption  permits  the  gas  dynamical  equations  to 
be  uncoupled  from  and  solved  independently  of  Ampere’s  law. 

The  assumption  has  been  used  extensively  in  magnetogasdynamic 
studies  and  is  equivalent  to  requiring  that  the  magnetic 
Reynolds  number,  R0  =  op.eUL,  be  small  compared  to  unity. 

With  these  assumptions  the  equations  of  motion  reduce 
to  the  following  form: 


Continuity  Equation: 


—  (pUj)  =  0  (12) 

dr 


r  momentum  equation: 
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pu. 


dr 


-  ^  d2U  A 

-  2E  +  nri  -  r  ■  4 


au 


dr 


+  • 


'3  Sr2  3r  dr 


3r 


1  S"Ur 

U  +  —  - -] 

2  1  ?  ae2  J 


+ 


<i 


—  U  )  +  — - -  lii  -  j  Bfl 

3r  r  dr  r2  d©  d©  Z  y 


(13) 


0  Momentum  Equation: 


0  = 


,  >  ,  d2U  ,,  dU 

I  ip  4-  li  (1 _ £  +  1 _ £) 


r  d0  3  r  drdS  r2  d© 


n-.au  U  >  ,  dU  N 

_  2  (I  _r  _  2  _r )  dji  +  I  dfi 

3  r  dr  r2  d9  r  59  dr 


(14) 


Energy  Equation: 


PUT 


dh  _ 

U  &  - 

dr 

r  dr 

r  au 

+  If 

|  2  ( — 

L  dr 

2 

- (X 

3 

Law: 

Jz 

dr 


2  Ur  2 
)  +  2<-£)  +  (■ 


r2  d©  d© 
dU. 


-) 


1  __r 

r  as 


au. 

f 

dr 


U  2 
-  +  — )  +  — 


J 

2 

a 


Jz  =  *<Ez  +  Ur8©) 


(15) 


(16) 


The  fluid  medium  is  assumed  to  behave  as  a  perfect 
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gas,  both  thermally  and  calorically  so  that: 


p  =  pRT 

(17) 

h  =  CpT 

(18) 

where  Cp  is  a  constant. 

To  complete  the  specification  of  the  problem,  boundary 
conditions  must  be  imposed  on  the  independent  variables 
Ur,  T,  and  p.  Along  the  channel  centerline,  0=0,  a  con¬ 
dition  of  symmetry  is  imposed. 


dUr(r,0) 

50 


=  0 


5T(r,Q)  =  Q 
50 


Zp(r,0)  =  0 

50 


(19) 


At  the  wall,  the  no  slip  condition  is  applied  and  the 
pressure  and  temperature  variations  are  specified. 


Ur(r,a)  =  0  ,  T(r,a)  =  Tw(r)  ,  p(r,a)  =  pw(r)  (20) 


Rather  than  specifying  the  additional  boundary  con¬ 
ditions  necessary  at  a  particular  value  of  r,  a  separable 
form  of  solution  is  assumed.  It  is  then  assumed  that  the 
dependent  variables,  transport  properties,  and  applied 
fields  can  be  expressed  in  the  following  form: 


Ur(r,6)  =  U(0)ra 
p(r,0)  =  p(9)rb 


(Equation  continued) 
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P (r , 0)  =  p (0)rC 
T(r,9)  =  T(9)rd 
Ez(r)  =  -  Ere 
Bg(r)  =  Brf 
p,(r,  0)  =  |I(0)rg 
o(r ,0)  =  o(0)rh 
k:  (r ,  6)  =  ic(0)r^ 


where  the  exponents  a,  b,  c,  etc.  are  constants  to  be 

determined.  Substitution  of  Eqs.  (21)  into  Eqs .  (12) 

through  (15)  and  eliminating  J  and  h  gives: 

z 


(a  +  b  +  l)pUra+b  =  0 


~7t2  sa+b-i  -  c-i  4  .  - 

apU^r  =  -  cpr  +  —  p,[a(a  -  1)U 

3 


+  au  -  U]ra+g-2  +  flu"ra+g-2  +  (i  aO  -  ^  B)g[Ira+g-2 

3  3 


+  U'p'ra+g  2  +  cEBre+f+h  -  c?UB2ra+2f+h 


0=-p'rCl+^  (au'  +  7u')ra+g-2 
3 


-  -  (aU  -  2U)d'ra+g"2  +  gU ,pra+g-2 
3 


(21) 


(22) 


(23) 


(24) 
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dpU<£Tra+b+d  1  =  cUpra+c_1  +  a(E  -  UB)2rh+2e 


+  d(d  +  j ) /cTrd+^-2  +  (kT"  +  /c'f')P+d  2 


+  \L  [2a2U2  +  2U2  +  <u' )2]r2a+g~2 


-  [L  (a  +  l)2U2r2a+S~2 


(25) 


prc  =  pRTrb+d 


(26) 


where  the  prime  refers  to  differentiation  with  respect 
to  6. 


If  U,  T,  etc,  are  to  be  functions  of  0  only,  as 
assumed,  then  Eqs .  (22)  through  (26)  may  not  involve  the 
variable  r.  Therefore  the  following  equations  must  be 
satisfied  by  the  exponents: 


a  +  b  +  1  -  0 


2a  +  b-  l  =  c-  l  =  a  +  g-  2 


=  e  +  f+  h  =  a  +  2f  +  h 


a  +  b  +  d-  1  =  a  +  c  -  1  =  h  +  2e 


=d+j  -2=2a+g-2 


c  -  b  +  d 


(27) 
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Eqs .  (27)  represent  only  seven  independent  equations 
for  the  nine  exponents  which  can  be  solved  in  terms  of  the 
electric  and  magnetic  field  exponents  e  and  f. 


a  =  e  -  f 

Ur  =  U (0)re-f 

b  =  f  -  e 

-  1 

p  = 

c  -  e  -  f 

-  1 

p  =  pie)!'0'1*1 

d  =  2  (e  - 

f) 

T  =  T(9)r£^e“f^ 

o 

ii 

bn 

d  =  il(0) 

h  =  -  2  (f 

+  1) 

a  =  5(9)r-2<f+l) 

C_l. 

II 

O 

K.  =  K  (9 ) 

(28) 


Since  the  fluid  transport  properties  cannot  vary  with 
r,  they  are  assumed  independent  of  8  as  well,  i.e.  p  = 
constant  and  k  -  constant  will  be  assumed. 

From  Eqs.  (28)  it  is  clear  that  the  radial  variations 
of  the  quantities  Ur,  p,  p,  and  T  are  determined  by  electric 
and  magnetic  field  distributions.  These  distributions  will 
be  chosen  so  as  to  be  consistent  with  Maxwell's  electro¬ 
magnetic  equations. 

The  r  and  8  components  of  Faraday's  law  for  steady 
state  are: 


c»(rEfl) 

- Z_  =  o 

90  9z 


9E  9E 

0  (29) 

9  z  9r 


Then  since  9/9z  =  0  by  assumption 


9r  99 


0 
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or 


E  =  constant,  e  =  0 
z 


(30) 


Such  an  electric  field  could  be  produced  by  using  parallel 
equipotential  electrode  surfaces  for  the  channel  sidewalls. 

A  solution  to  Ampere's  law,  Eq.  (6)  is 

B  =  Bh  +  b  (31) 


where 


V  x  Br  =  0 


-> 

V  x  b  =  |ieJ 


It  is  clear  that  BR,  the  solution  to  the  homogeneous 

equation,  is  the  applie<|  field  Bg  and  that  b  is  the  induced 

field.  By  assumption  |b|  «  Bg.  Thus  the  applied  magnetic 

field  Bg  should  satisfy  the  z  component  of  Ampere’s  law 

with  J  =0. 
z 


or 


S(rB0) 

Sr 


(32) 


Such  a  magnetic  field  could  be  produced  by  a  current  carrying 
wire  coincident  with  Z  axis  and  could  doubtless  be  closely 
approximated  by  a  magnet  with  diverging  pole  faces. 

It  should  be  pointed  out  that  while  Eqs.  (30)  and  (32) 
define  the  only  electric  and  magnetic  field  distributions 
rigorously  consistent  with  the  geometry  and  assumptions  of 
this  model,  in  an  actual  three-dimensional  channel  other 
distributions  are  not  only  possible  but  may  be  desirable. 
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By  using  finely  segmented,  individually  powered  electrodes 
almost  any  reasonable  electric  field  distribution  can  be 
closely  approximated.  Although  the  magnetic  field  is  not 
so  easily  modified,  this  distribution  too  can  be  tailored 
to  some  extent. 

Introducing  the  values  of  e  and  f  into  Eqs.  (28) 
gives : 


a  =  1  Up  =  U(0)r 

b  =  -  2  p  =  p(0)r 


c  =  0  p  =  p(9) 

l 

d  =  2  T  =  T(S)r2 


h  =  0 


a  -  a(0)  =  constant  (33) 


Since  the  electrical  conductivity  cannot  vary  with  r,  it 
will  be  taken  to  be  a  constant. 

Eqs.  (33)  show  that  the  gas  is  accelerated  down  the 
channel,  the  pressure  and  Mach  number  being  constant  along 
radial  lines.  It  is  not  surprising  that  these  radial 
variations  are  of  the  same  form  as  found  by  Hains  (Ref.  8) 
and  Podolsky  and  Borman  (Ref.  9)  on  the  basis  of  inviscid 
models. 

Substituting  the  values  of  the  exponents  into  Eqs. 
(22)  through  (28)  gives: 


pU2  =  |iu"  +  a(E  -  UB)B 


(34) 


P 


(35) 


2pUC  T  =  c(E  -  UB)2  +  kT"  +  4k: T 
P 

+  [l[  <U')2  4-  -  U2] 

3 


(36) 
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P  =  pRT 


(37) 


Eq.  (22)  is  identically  satisfied. 

In  summary  the  additional  restrictions  that  have  been 
forced  by  the  separation  of  variables  are: 

1.  The  variations  of  velocity,  density,  temperature 
and  pressure  down  the  channel  have  been  specified. 

2.  It  has  been  shown  that  the  viscosity,  thermal 
conductivity,  and  electrical  conductivity  can¬ 
not  change  with  r.  These  properties  have  been 
taken  to  be  constant. 

3.  The  variation  of  temperature  and  the  constancy 

of  pressure  have  been  specified  along  the  channel 
wall. 


T(r,cc)  =  Tw(r)  =  T^r2 


p(r,a)  =  pw(r)  =  pw 


Eq.  (35),  the  6  momentum  equation,  can  be  integrated 
immediately  to  give 


p  =  £  [LV  +  p  (38) 

3  w 

Using  Eq.  (38)  and  the  equation  of  state  given  by  Eq.  (22) 
the  density  p  can  be  eliminated  from  Eqs.  (34)  and  (36) 
giving: 


—  (-  M-U  +  p  )  =  (lU"  +  o(E  -  UB)B  (39) 

RT  3 


R 


0 (E  -  UB)2  +  icT"  +  4kT 

(Equation  continued) 


20 


AE  DC-TR-67-80 


+  ^[(U')2  +  -  (1  -  ^£)U2] 
3  R 


(40) 


Eqs.  (39)  and  £40)  are  two  coupled,  non-linear  differential 
equations  for  U  and  T.  The  boundary  conditions  for  these 
two  quantities  are: 


U'(0)  =  0  U (a)  =  0 

T'(0)  =  0  T(a)  =  Tw  (41) 


Eqs.  (39)  and  (40)  can  be  put  into  dimensionless  form  by 
defining  new  variables 


f(0) 


g(9) 


(42) 


where  rQ  is  the  radial  distance  to  the  channel  "entrance", 
Uc  and  T„  are  characteristic  values  of  velocity  and  tempera 
ture  to  be  specified.  Also  define 


B 

o 


(43) 


where  B0  is  the  magnetic  induction  at  rG.  Define  the 
dimensionless  parameters 


Ha 


<-M  B 
M- 


r 


o 


o 


Re 


p_U  r 
rc  c  o 


M- 


(Hartmann  number) 


(Reynolds  number) 


where 


c 


(44) 
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9  ~ 


U  B 
c  o 


uc 

M  =  - t  (Mach  number) 

(7RTc)2 


C  \l 

Pr  =  — (Prandtl  number) 


Then  Eqs .  (39)  and  (40)  can  be  rewritten  in  terms  of  the 
dimensionless  variables  and  parameters  as 


=  I  7m2  -3  +  Re  -S  -  (Ha)2 ($  -  f) 
d02  3  g  g 


2 

+  4g  =  2RePrf  -  (7  -  l)M2Pr  (Ha) 2  (<3>  -  f)2 

d02 


-  (7  -  l)M2Pr(— )2  +  —  (37  +  l)M2Prf2 
d  9  3 


with  boundary  conditions 

f'(0)  =  0 

g'(0)  =  0 


Eqs.  (45)  and  (46)  with  boundary  conditions  given  by 
Eq.  (47)  represent  the  analytical  model  for  the  present 
investigation. 


f (a)  =  0 


g(a)  = 


T  r 
w  o 


(45) 


(46) 


(47) 
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IV 

TECHNIQUE  OF  SOLUTION 


4.1  SOLUTION  DIFFICULTIES  OF  THE  DIRECT  BOUNDARY  VALUE  PROBLEM 


Eqs.  (45)  through  (47)  appear  to  define  a  straight 
forward  but  complicated  non-linear  boundary  value  problem. 

The  following  method  of  solution  would  seem  applicable. 

1.  Assign  values  for  the  divergence  angle  a 
and  the  wall  temperature  g(a). 

2.  Specify  values  for  the  six  dimensionless 
parameters  7,  Pr,  Ha,  Re,  M,  $. 

3.  Guess  initial  values  for  the  dimensionless 
velocity  and  temperature  f(0)  and  g(0). 

4.  Numerically  integrate  Eqs.  (45)  and  (46) 
using  the  Runge-Kutta  method  from  0=0 
to  9  -  a  hoping  to  satisfy  the  boundary 
conditions  at  the  wall. 

5.  Adjust  the  values  of  f(0)  and  g(0)  as 
necessary  and  integrate  the  equations 
again.  Repeat  until  the  boundary  con¬ 
ditions  at  a  are  met. 

6.  Assign  new  values  for  the  parameters  7, 

Pr,  Ha,  Re,  M,  0  and  repeat  steps  3  through 

5. 

The  iterative  procedure  outlined  above  would  obviously 
involve  considerable  computation  time  with  no  a  priori 
guarantee  of  convergence.  There  is  also,  however,  another 
somewhat  more  subtle  difficulty  with  this  method. 

In  order  that  the  dimensionless  variables  and  parameters 
be  well  defined,  the  characteristic  velocity  Uc  and  the 
characteristic  temperature  Tc  must  be  specified.  A  convenient 
value  for  the  characteristic  temperature  is  the  wall  temper¬ 
ature  at  rQ,  Tc  =  Twrp2,  so  that  the  temperature  boundary 
condition  becomes  g (a)  =  1.  For  the  characteristic  velocity, 
a  logical  choice  would  be  the  average  velocity  at  rQ,  defined 
as 

r  a 

U  —  — ^  f  U(0)d6  (48) 
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Introducing 


f(0) 


U(9)r. 


U 


gives 


a 

/  f(0)d 8  =  2a  (49) 

-a 


Def iningvr)  =  8/a  and  using  the  fact  that  f(G)  is  an  even 
function  gives 

01  1 

/f(n)dr|  =  1  (50) 

0 


an  extra  condition  of  f  that  generally  would  not  be  satisfied. 
For  problems  that  can  be  solved  explicitly  in  closed  form, 
this  extra  condition  usually  defines  a  relationship  between 
the  parameters.  For  example,  in  ordinary  hydrodynamic  flow 
between  parallel  plates  a  condition  similar  to  Eq.  (50)  gives 
a  relation  between  the  mean  velocity  and  the  pressure  gradient. 

Thus  in  the  present  problem  one  would  expect  some 
relationship  among  the  parameters  y ,  Pr,  Ha,  Re,  M,  and 
that  must  be  satisfied  before  a  solution  can  be  obtained. 

This  conclusion  can  be  illustrated  by  considering  the 
corresponding  one-dimensional,  inviscid  problem. 

The  governing  equations  for  the  inviscid  problem  can 
be  obtained  from  Eqs.  (12)  through  (18)  by  setting  p.  =  k  -  0. 
The  resulting  equations  are 


3(pU  r) 

- £ —  e  o 

dr 


(51) 


-  ^  +  a(E  -  U  B)B 
dr  r 


(52) 
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pU  C  — 
P  r  P  dr 


=  U  ^ 


Sr 


4-  O  (E  — 


U  B)2 
r 


(53) 


p  =  pRT 


(54) 


Define  new  dimensionless  variables  as 


where  the  subscript  zero  denotes  the  value  of  the  quantity 
at  rQ.  The  applied  magnetic  field  distribution  is  the 
same  as  previously  derived  in  Eq.  (32),  i.e.  B  =  B0/r* . 

In  terms  of  the  dimensionless  variables  the  equations 
become 


p*U* 


dU* 

dr* 


ir 


V*  - 


(56) 


p*U*r*  -  1 


(57) 


p*U*  dT* 
(7-1 )M2  dr* 


irU*  &L. 
dr* 


+  R  (<£>  - 
M 


2 

r 


(58) 


p*  =  p*T* 


(59) 
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where 


<t  = 


UB 
o  o 


*11 


or  B  2 
o  o 

P  U 
“  o  o 


7 r  = 


p  U  " 
^  o  o 


U 


M  = 


(7RT0)2 


(60) 


The  initial  conditions  on  the  dependent  variables  are 

U*(l)  =  1 
T*(l)  =  1 
p*(l)  =  1 

P*(l)  =  1  (61) 

Following  a  procedure  similar  to  that  used  for  the  two 
dimensional  problem  assume  a  solution  of  the  form 

U*  =  Cx(r*)a 
p*  =  C2(r*)b 
p*  =  C  (r*)c 

3 

T*  =  C  (r*)d  (62) 
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which  leads  to 


u*  M, 


1  +  Si 


1  +  R.i  _ 

p*  =  - 1  (r*)"2 

R  4> 

Sr 


* 

P  = 


(7-1  )M2 

2<J  +  V 


T*  ,  (^-i)Ma¥a 


2(1  +  V 


(63) 


(64) 


(65) 


(66) 


This  solution  will  satisfy  the  initial  conditions  if  and 
only  if, 


y 

1  +  Si 


1 


(7  -  1)M2R1I02 

2(1  +  V2  1 


or 


_  1  +  Si 

Si 


(67) 


M2 


2Si 

7-1 


(68) 


If  Eqs.  (67)  and  (68)  are  satisfied  the  solution  is 
given  by 
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U*  =  r* 

P*  =  (r*)~2 
P*  =  1 

T*  =  (r*)£  (69) 


If  Eqs.  (67)  and  (68)  are  not  satisfied  the  solution 
is  entirely  different,  but  according  to  Podolsky  and  Borman 
(Ref.  9)  approaches  the  values  given  by  Eqs.  (63)  through 
(66)  for  large  values  of  r*.  The  solution  given  above  is 
a  special  case  of  the  problem  investigated  by  these  authors 
although  it  was  not  specifically  pointed  out  by  them. 

In  view  of  the  above  considerations  one  is  forced  to 
conclude  that  the  direct  boundary  value  problem  is  indeed 
overspecified  and  that  a  solution  cannot  be  obtained  for 
an  arbitrary  choice  of  the  parameters.  This  overspecification 
would  make  solution  by  the  iterative  procedure  outlined 
previously  extremely  difficult.  Fortunately,  there  is  an 
alternate  approach  to  the  problem  that  alleviates  this 
difficulty. 

4.2  THE  INITIAL  VALUE  PROBLEM 

An  approach  which  will  circumvent  the  difficulty  of 
overspecification  of  the  problem  mentioned  above  is  to 
leave  the  wall  divergence  angle  as  a  dependent  variable  of 
the  problem.  The  problem  also  becomes  an  initial  value 
problem  in  contrast  to  the  two  point  boundary  value  problem 
specified  by  Eqs.  (47). 

For  this  approach  the  characteristic  velocity  and 
characteristic  temperature  of  the  problem  are  defined  to 
be  values  of  velocity  and  temperature  at  r  =  r  and  0=0, 
i.e.  ° 


U  =  U  =  U(r  ,0) 
c  o  o 

T  =  T  =  T(r  ,0)  (70) 

coo 


With  this  choice  the  dimensionless  equations  and  boundary 
conditions  become 
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flS  J!  Q  _  f-3  fS 

=  -  7M2  —  +  Re  - - (Ha)2 ($  -  f)  (71) 

d02  3  g  g 


^  +  4g  =  2RePrf  -  (7  -  l)M2Pr (Ha)2($  -  f)2 
d02 


(7  -  l)M2Pr  (— )  +  -  O7  +  1 )M2Prf 2 

d@  3 


(72) 


f(0)  =  1 
g(0)  =  1 


f'(0)  =  0 


g  (0)  =  0 


(73) 


Eqs.  (71),  (72)  and  (73)  constitute  an  initial  value 
problem  that  can  be  integrated  directly  by  numerical  methods. 
The  divergence  angle  and  the  wall  temperature  are  then 
determined  by  the  integration.  The  no  slip  condition 
f(a)  —  0  is  used  to  specify  the  wall  divergence  angle,  a, 
and  the  corresponding  value  g(a)  gives  the  ratio  of  wall 
to  centerline  temperatures. 

Thus,  in  approaching  the  problem  by  the  method  des¬ 
cribed  above,  two  new  parameters  have  been  introduced, 
giving  a  total  of  eight,  namely  7,  Pr,  Ha,  Re,  M,  a, 
g(a).  In  order  to  keep  the  number  of  computations  required 
within  reason,  it  was  decided  to  fix  and  hold  constant  the 
ratio  of  specific  heats,  7,  and  the  Prandtl  number  at  values 
typical  for  diatomic  gases.  Thus  the  values  7  =  1.4  and 
Pr  =  1  are  used  in  all  numerical  calculations. 

Eqs.  (71)  and  (72)  with  initial  conditions  given  by 
Eqs.  (73)  were  programmed  for  solution  using  a  standard 
fourth  order  Runge-Kutta  scheme. 

Consider  now  the  interrelation  among  the  remaining 
parameters  Ha,  <t,  Re,  M,  a.,  g(a).  It  is  to  be  expected 
that  a  change  in  any  one  of  the  parameters  Ha,  <±>,  Re,  or  M 
will  cause  a  change  in  both  a  and  g(a).  Extensive  computer 
calculations  reveal  that  this  is  true.  More  generally  it 
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may  be  said  that  only  three  of  the  six  parameters  may  be 
held  constant  at  once  or  that  a  change  in  any  one  of  these 
must  result  in  a  change  in  two  of  the  remaining  five. 

Figs.  3  and  4  illustrate  this  behavior  showing  the  relation 
between  Ha,  <i>,  a,  and  g(a)  for  Re  and  M  fixed.  Lines  of 
constant  a  have  been  cross  plotted  on  Fig.  4. 

Thus  it  is  not  possible  in  this  problem  to  directly 
assess  the  effect  of  the  change  of  a  single  parameter. 

Other  solutions  for  viscous,  compressible,  channel  flows 
both  with  and  without  electromagnetic  effects  exhibit  a 
similar  behavior  (Ref.  2,  19,  20). 

In  view  of  these  considerations,  solutions  will  be 
compared  on  the  following  basis 

1.  The  electromagnetic  parameters  Ha  and  $  will 
be  taken  as  the  independent  parameters. 

2.  The  values  a  and  g(a)  will  be  considered 
as  fixed. 

3.  The  gas  dynamic  parameters  Re  and  M  will  be 
considered  dependent  parameters. 

This  technique  will  permit  a  comparison  of  velocity  and 
temperature  profiles  for  flows  with  Ha  =  Ha)-,^  and  Ha  =  Ha)2 
at  the  same  values  of  <t>,  a,  and  g  (a)  .  Corresponding  to 
the  change  in  Ha  there  will  be  a  change  in  both  Re  and  M. 

It  is  felt  that  this  method  will  yield  the  most  information 
about  the  effect  of  a  change  in  one  of  the  electromagnetic 
parameters  on  the  flow. 

From  the  numerical  calculations,  graphs  were  constructed 
showing  a  as  a  function  of  Re  and  M  and  g(a)  as  a  function 
of  Re  and  M  for  fixed  values  of  Ha  and  $. 

Range  of  Parameters 

It  was  quickly  discovered  through  numerical  computations 
that  there  are  well  defined  limits  on  the  values  that  some 
of  the  parameters  can  be  assigned.  For  example,  reference 
to  Fig.  26  shows  that  there  is  a  minimum  value  of  the 
divergence  angle  which  occurs  at  Re  =  0.  Moreover,  it  is 
clear  that  (da/dM)Re  _  g  =  0.  From  these  observations  an 
expression  can  be  obtained  for  the  minimum  divergence  angle 
as  a  function  of  Ha  and  $>.  Putting  Re  =  0  and  M  =  0  in 
Eq.  (45)  gives 

—  +  (Ha)*($  -  f)  =  0  (74) 

d02 
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with  initial  conditions 


f (0)  =  1  f  '  (0)  =  0 


(75) 


which  yields 


f  =  4>  4-  (1  -  ®)cosh  (Ha9) 


(76) 


Using  the  condition  f(a)  =  0  gives 


a 


min 


—  cosh 
Ha 


(77) 


For  Ha  =  100,  $  =  1.6,  Eq.  (77)  gives  am^n  -  .0164  rad  =  .94° 
which  agrees  with  Fig.  28.  Also  the  general  behavior 
predicted  by  Eq.  (77)  is  the  same  as  shown  by  Fig.  3  for 
non-zero  values  of  Re  and  M. 

Note  that  as  4>  approaches  unity  amin  approaches  infinity 
and  that  for  $  less  than  one,  no  real  values  of  amin  exist. 
This  would  suggest  that  meaningful  solutions  to  Eqs.  (45) 
and  (46)  are  only  possible  for  $  >  1 .  This  conclusion  has 
been  verified  by  numerical  calculations. 

The  quantity  =  E0/UoB0  is  an  electromagnetic  loading 
parameter  and  0  <  4>  <  1  corresponds  to  a  generator  mode 
of  operation  of  the  channel  where  energy  is  extracted  from 
the  fluid.  The  accelerator  mode  of  operation  corresponds 
to  $  >  1.  Since  the  solution  form  assumed  in  Eqs.  (33) 
clearly  indicates  that  energy  is  being  added  to  the  fluid 
by  the  applied  fields,  it  is  not  surprising  that  solutions 
cannot  be  obtained  for  $  <  1. 

Having  obtained  an  expression  for  the  minimum  value  of 
the  divergence  angle,  it  is  logical  to  inquire  whether  or 
not  this  parameter  has  an  upper  limiting  value.  Rewriting 
the  differential  equations  for  f(8)  and  g(0)  and  the  initial 
conditions  on  these  variables  gives 


—  =  ®  yM2  —  4-  Re  —  -  (Ha)2($  -  f)  (78) 

d  0Z  3  g  g 
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+  4g  =  2RePrf  -  (7  -  1  )M2Pr  (Ha)2  ($>  -  f)2 

d02 


-  (7  -  l)M2Pr(— )2  +1  (37  +  l)M2Prf 2  (79) 

d0  3 

f(0)  =  1  f'(0)  -  0 

g(0)  *=  1  g '  (0)  =  0  (80) 

Then  at  9  =  0 

f  "(0)  =  -  7M2  +  Re  -  (Ha)2(3>  -  1)  (81) 

3 


g"(0)  =  2RePr  -  4  -  (7  -  1  )M2Pr  (Ha) 2  ($>  -  l)2 


+  -  (37  +  l)M2Pr 
3 


(82) 


If  Re,  M,  Ha,  and  <5  are  chosen  such  that  f"(0)  =  0  and 
g"(0)  =  0,  it  is  seen  by  continued  differentiation  of 
Eqs.  (78)  and  (79)  and  setting  0=0  that 


f n(0 )  =  0  gn(0)  =  0  n-3 , 4 . . .  (83) 


Thus 

f(0)  =  1  g(0)  =  1 


for  all  values  of  0  and  a  =  00.  Obviously  then  there  is  no 
upper  limit  on  the  values  that  a  can  attain.  Nonetheless, 
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a  very  interesting  result  has  been  obtained.  Considering 
Ha  and  $  as  being  fixed,  the  equations  corresponding  to 
f"(0)  =  0  and  g"(0)  =  0,  namely 


-  7M2  +  Ee  -  (Ha)2($  -  1)  =  0  (84) 

3 


2RePr  -  4  -  (7  -  l)M2Pr  (Ha)2(<£  -  l)2 


+  -  O7  +  l)M2Pr  =  0 
3 


(85) 


define  a  limiting  point  in  the  Re,  M  plane. 

It  is  expected  that  for  fixed  values  of  Ha  and  as 
values  of  Re  and  M  are  chosen  approaching  the  limiting 
values,  defined  by  Eqs.  (84)  and  (85),  larger  values  of  a 
will  be  obtained. 

It  is  particularly  interesting  to  note  the  close 
relationship  between  Eqs.  (84)  and  (85)  and  Eqs.  (67)  and 
(68),  the  conditions  that  must  be  satisfied  in  order  that 
the  inviscid  solution  be  given  by  Eqs.  (69).  If  the  terms 
8  4 

tt  7M  ,  4,  and  -r  (37  +  l)M2Pr  are  ignored,  Eqs.  (84)  and  (85) 
become 


Re  -  (Ha) 2  (<J>  -  1)  =  0 


(86) 


2RePr  -  (7  -1 )M2Pr (Ha) 2 (<f>  -  l)2  =  0  (87) 

But 

(Ha)2  =  —  B  2r  2 
M-  °  ° 

so  that  Eqs.  (86)  and  (87)  become 


or  B  2  p  U  r 
=  (  0  0  )(  °  0  < 


P  U 

o 


M- 


:)  =  HMRe 
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1  +  R*. 

<i>  =  - —  (88) 

% 


M 


2 


2rm 

7-1 


(89) 


Eqs.  (88)  and  (89)  are  identical  to  Eqs.  (67)  and  (68). 

The  terms  ^  7M2,  4,  and  —  (37  +  l)M2Pr  arise  due  to 
o  3 

the  cross  channel  pressure  variation,  axial  heat  conduction, 
and  viscous  dissipation  due  to  axial  velocity  gradients. 

If  the  problem  had  been  formulated  not  including  these 
effects,  i.e.  using  boundary  layer  type  approximations  to 
simplify  the  equations,  then  the  inviscid  solution  given 
by  Eqs.  (69)  would  represent  truly  the  limiting  solution 
for  the  viscous  problem.  With  the  aforementioned  effects 
excluded  the  dimensionless  equations  for  f(0)  and  g(0)  are: 


=  Re  —  -  (Ha)2 (0  -  f)  (90) 

d  9s  g 


=  2RePrf  -  (7  -  l)M2Pr (Ha) 2 (<J>  -  f)2 

d02 


-  (7  -  1  )M2Pr  (— )2 
dO 


with  initial  conditions 


f(0)  =  1 

g(0)  -  1 


f'(0)  =  0 
g'  (0)  =  0 


(91) 


(92) 


Setting  0=0  gives 
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—  f"(0)  =  1  -  IL.($  -  1)  (93) 

Re 


— g"(0)  =  2  -  {7  -  1)M2R ;  ($  -  l)2  (94) 

RePr 


so  that  for 


a  1  +  *M 


M2 


(95) 


fn(0)  =  gn(0)  =0  n  =  2,3,... 


(96) 


f<0)  =  g(©)  =  1 


(97) 


for  all  9  and  for  all  values  of  Re  and  Pr.  Then  the 
solution  is 


U  r(r,0)  =  U0(— ) 
r 


(98) 


T(r,9)  =  T  (— )2 
r 


(99) 


This  is,  of  course,  the  same  solution  obtained  for  the 
inviscid  problem  and  is  valid  only  when  Eqs .  (67)  and  (68) 
are  satisfied. 

Returning  to  a  consideration  of  the  limits  on  the 
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various  parameters,  it  was  decided  to  concentrate  on  values 
of  g(a)  between  0.1  and  0.5.  It  is  unlikely  that  the  ratio 
of  wall  temperature  to  centerline  temperature  would  be 
outside  these  limits  in  a  practical  device.  Fig.  5  shows 
the  region  in  the  Re,  M  plane  where  acceptable  solutions 
can  be  obtained  for  Ha  =  200  and  O  =  1.6.  As  was  anticipated, 
as  values  of  Re  and  M  are  chosen,  in  this  region,  approaching 
the  limit  point,  larger  values  of  divergence  angle  are  ob¬ 
tained.  Note  that  the  region  of  acceptable  solutions  becomes 
increasingly  narrow  as  the  limiting  point  is  approached. 

If  values  of  Re  and  M  are  used  that  lie  to  the  left 
of  the  region  shown  in  Fig.  5,  negative  values  of  g(a)  are 
obtained.  Such  a  result  is  physically  meaningless.  If 
values  are  used  to  the  right  of  the  region  shown,  large 
values  of  g(a)  are  obtained  and  ultimately  the  dimensionless 
velocity,  f<6),  grows  without  limit. 

Finally,  there  are  practical  upper  limits  on  the  values 
of  Ha,  and  a.  As  is  shown  in  Fig.  3  and  by  Eq.  (77),# 
a  tends  to  decrease  with  increasing  values  of  Ha  and  $. 

Thus,  for  large  values  of  these  parameters,  it  is  necessary 
to  closely  approach  the  limiting  values  of  Re  and  M  in  order 
to  hold  a  constant.  But,  as  was  pointed  out,  the  region  of 
acceptable  solutions  becomes  very  narrow  in  the  neighborhood 
of  the  limiting  point.  Thus,  it  is  very  difficult  to  obtain 
solutions  for  large  values  of  Ha,  $,  and  a  simultaneously. 

In  view  of  these  considerations  it  was  decided  to  con¬ 
centrate  on  values  of  a  between  1.0°  and  1.6°,  to  use  values 
of  Ha  =■  100,  150,  200,  and  to  use  values  of  4>  =  1.2,  1.6, 

2.0.  For  these  selections  it  is  possible  to  carry  out  the 
proposed  solution  and  comparison  scheme  with  a  reasonable 
amount  of  numerical  calculations. 

In  summary,  it  has  been  ascertained  that  there  are 
well  defined  limitations  on  the  range  of  each  of  the  para¬ 
meters  involved  in  the  present  study,  Ha,  Re,  M,  a,  and 
g(a).  These  restrictions  are  either  mathematical,  physical, 
or  practical,  depending  on  the  parameter  but  in  all  cases 
they  are  quite  real  and  meaningful  solutions  can  be  obtained 
only  within  these  limitations.  Table  1  gives  the  numerical 
values  of  amj_n,  ^eiimit  ^or  ^he  values  of  Ha  and 

used  in  this  study. 
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V 

DISCUSSION  OF  RESULTS 


The  results  of  this  portion  of  the  investigation  are 
presented  graphically  in  Figs.  6  through  18.  Figs.  6  through 

14  give  velocity  and  temperature  profiles  for  the  nine  com¬ 
binations  of  Ha  and  <t>  considered  in  the  study.  The  dimen¬ 
sionless  variables  f  and  g  are  plotted  against  the  normalized 
space  coordinate  rj  =  0/a4  For  each  of  these  graphs  the 
divergence  angle  is  approximately  1.5°  and  the  ratio  of  wall 
temperature  to  centerline  temperature  is  approximately  0.2. 

The  values  of  Re  and  M  required  to  maintain  these  values 

a  and  g(a)  were  obtained  from  Figs.  26  through  43.  Figs. 

15  and  16  show  the  variation  of  viscous  drag  and  wall  heat 
transfer  with  Ha  and  4>.  Figs.  17  and  18  are  velocity  and 
temperature  profiles  for  a  =  1.25°  and  g(a)  =  0.2.  These 
figures  show  the  effect  of  a  change  in  the  divergence  angle 
and  will  be  used  for  comparison  with  the  results  for  variable 
properties . 

The  velocity  profiles  can  also  be  interpreted  to  give 
the  cross  channel  pressure  variation.  Writing  Eq.  (38)  in 
terms  of  dimensionless  variables .and  parameters  and  letting 


P*(S) 


P<9)  -  P„ 


gives 

p*(G)  =  ®  ^  f(0) 
3  Re 


(100) 


(101) 


Table  II  summarizes  the  results  giving  the  values  of 
Re,  M,  f'(ct),  g'(a),  and  p*(0)  in  terms  of  Ha  and  $  for 
a  =  1.5°  and  g(a)  =  0.2. 

It  must  be  pointed  out  again  that  the  changes  in  the 
velocity  and  temperature  profiles  exhibited  in  Figs.  6 
through  14  cannot  be  solely  attributed  to  changes  in  the 
electromagnetic  parameters  Ha  and  <P.  As  has  been  discussed, 
associated  with  a  change  in  one  or  both  of  these  parameters, 
there  is  a  change  in  the  gas  dynamic  parameters  Re  and  M. 
This  interdependence  of  electromagnetic  and  gas  dynamic 
parameters  will  be  shown  in  the  following  discussion  to  be 
an  important  feature  of  the  problem. 
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Perhaps  the  most  salient  feature  of  the  results  is 
the  prominent  bulge  in  the  temperature  profiles.  Physically 
this  is  due,  of  course,  to  the  concentration  of  joule  heating 
and  viscous  heating  in  the  vicinity  of  the  cool  wall.  In 
this  region  the  reversed  electromotive  force  due  to  the 
fluid  motion  is  small  and  a  large  current  density  results. 
This  effect,  however,  is  somewhat  more  pronounced  than  had 
been  anticipated. 

The  trend  exhibited  by  the  temperature  bulge  as  <I> 
changes  is  very  interesting.  If  gmax  denotes  the  maximum 
value  of  g(t|),  for  certain  values  of  Ha,  gmax  increases 
monotonically  with  <&.  At  a  larger  value  of  Ha,  gmax  in- 
creases  to  a  maximum  value  and  then  decreases  with  increasing 
$,  and  for  a  still  larger  value  of  Ha,  gmax  may  monotonically 
decrease  with  <E>.  In  particular  for  Ha  -  100,  gmax  increases 
monotonically  with  0,  for  Ha  =  150,  gmax  increases  to  a 
maximum  and  then  decreases  with  $,  and  for  Ha  =  200,  gmax 
decreases  monotonically  with  <5.  This  behavior  is  also 
reflected  in  Fig.  16  which  shows  the  temperature  gradient 
at  the  wall  as  a  function  of  $  for  the  above  three  values 
of  Ha. 

In  order  to  explain  this  phenomenon,  the  behavior  of 
the  temperature  profile  for  constant  values  of  Ha  and 
along  a  line  of  constant  g(a)  in  the  Re,  M  plane  must  be 
examined.  Now  for  g(a)  fixed  there  is  a  correlation  between 
the  bulge  in  the  temperature  profile  and  g"(0). 


g"(0)  =  2RePr  -  4  -  (7  -  l)M2Pr (Ha)2 (4>  -  l)2 


+  -  (37  +  l)M2Pr  (102) 

3 


From  Eq.  (102),  considering  the  dominant  first  and  third 
terms  on  the  right  side,  it  is  clear  that  g"(0)  increases 
with  increasing  Re  and  decreases  with  increasing  values  of 
M.  Thus,  referring  to  Fig.  5,  one  would  expect  upon  moving 
along  a  line  of  constant  g(ct)  toward  the  limit  point  that 
the  bulge  in  the  temperature  profile  would  first  increase, 
reach  a  maximum  and  then  decrease.  The  temperature  profile 
must  flatten  as  the  limit  point  is  approached. 

Next  consider  a  three-dimensional  O,  Re,  M  space  for 
fixed  Ha.  In  such  a  space  there  are  surfaces  of  constant 
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a  and  surfaces  of  constant  g(a).  Corresponding  to  the 
limiting  point  in  the  Re,  M  plane,  there  is  a  limiting 
curve  in  0,  Re,  M  space.  In  obtaining  solutions  for  Ha, 
a,  and  g(a)  fixed  and  for  various  values  of  0,  one  obtains 
solutions  corresponding  to  points  in  this  space  that  lie 
on  the  line  of  intersection  of  a  constant  g(a)  surface  and 
a  constant  a  surface.  Now  it  was  stated  in  Section  IV  that 
for  larger  values  of  $  it  is  necessary  to  closely  approach 
the  limiting  values  of  Re  and  M  in  order  to  hold  a  constant. 

In  other  words  as  the  coordinate  increases  the  line  of 
intersection  approaches  the  limiting  curve.  Thus  in  obtaining 
solutions  on  the  line  of  intersection  for  increasing  i.e. 
approaching  the  limiting  curve,  one  expects  a  behavior 
similar  to  that  observed  in  moving  along  a  line  of  constant 
g(a)  in  the  Re,  M  plane,  that  is  an  increase  in  the  temper¬ 
ature  bulge  followed  by  a  decrease. 

Mathematically,  then,  the  reason  for  the  aforementioned 
behavior  becomes  clear.  For  Ha  =  100  increasing  &  from  1.2 
to  2.0  along  the  line  of  intersection  of  the  surfaces  a  =  1.5° 
and  g(a)  =  0.2  occurs  on  a  section  of  this  line  removed 
far  enough  from  the  limiting  curve  such  that  the  temperature 
bulge  becomes  more  pronounced.  For  Ha  =  200  the  change  in 
<i>  occurs  on  a  section  of  the  corresponding  line  close  enough 
to  the  limiting  curve  so  that  the  temperature  bulge  becomes 
less  pronounced. 

While  the  explanation  given  above  is  perhaps  satisfactory 
from  a  mathematical  point  of  view  it  yields  little  insight 
into  the  physical  mechanisms  involved. 

The  dimensionless  parameter  $  can  be  interpreted  as  a 
measure  of  the  total  energy  input  divided  by  the  accelerating 
work . 


Total  energy  input/ (unit  vol)(unit  time)  =  J  ♦  ? 
Accelerating  work/ (unit  vol)(unit  time)  -  (J  x  B) 
Forming  the  ratio  of  these  quantities  gives 


-> 

U 


where 


J  •  E 


J  E 
z  z 


$ 


(J  x  £)  .  tJ  JzB0Ur  f(0) 


$  = 


BU 
o  o 


(103) 
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But 

->  ->  t2  ->  ->  -> 

J  »  E  =  — —  +  (J  x  B)  •  U 

a 

so  that  Eq,  (103)  can  be  written  as 

-  J!/0  -  +  1  (104) 

f(0)  (J  x  B)  •  U 


Eq.  (104)  indicates  that  for  larger  values  of  0  a  greater 
portion  of  the  total  energy  input  is  dissipated  in  Joule 
heating.  It  has  been  noted  that  this  dissipation  is  largest 
in  the  vicinity  of  the  wall  where  the  reversed  electromotive 
force  due  to  the  fluid  motion  approaches  zero.  Then  as 
expected,  for  low  values  of  Ha  this  increased  dissipation 
for  larger  values  of  $  is  exhibited  as  a  more  pronounced 
temperature  bulge. 

Evidently,  however,  another  mechanism  is  involved  that 
counteracts  this  effect  and  dominates  it  for  large  values 
of  the  Hartmann  number.  Examination  of  the  velocity  profiles 
shows  that  there  is,  as  expected,  a  flattening  of  these 
profiles  with  increasing  Ha  and  to  a  lesser  degree  with 
increasing  $.  It  is  suggested  that  this  flattening  of  the 
velocity  profiles  for  large  values  of  the  Hartmann  number 
tends  to  depress  the  temperature  bulge  by  distributing  the 
Joule  dissipation  more  evenly  over  the  channel. 

Obviously,  the  physical  explanation  given  above  is  not 
the  complete  story;  the  effects  of  viscous  dissipation  and 
heat  transfer  to  the  wall  certainly  play  an  important  role 
in  the  behavior  of  the  temperature  profiles.  It  is  felt, 
however,  that  the  dissipation  of  energy  in  Joule  heating 
is  the  most  important  mechanism  affecting  the  temperature 
distribution  in  the  channel. 

Another  interesting  feature  of  the  results  is  the  slight 
bulge  in  the  velocity  profiles  that  occurs  at  the  large  values 
of  Ha  and  <J>.  This  behavior  is  no  doubt  due  to  the  fluid 
density  variation  across  the  channel.  Table  II  shows  that 
the  pressure  change  across  the  channel  is  very  small  so 
that  the  density  is  very  nearly  inversely  proportional  to 
the  temperature.  Thus  regions  of  high  temperature  correspond 
to  regions  of  low  fluid  density.  Now  the  J  x  B  body  force 
is  a  volumetric  force  and  therefore  strongly  affects  the 
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fluid  in  regions  of  low  density.  It  is  concluded  then  that 
the  bulge  in  the  velocity  profiles  is  caused  by  the  density 
variation  across  the  channel  coupled  with  the  electromagnetic 
body  force.  The  fact  that  the  velocity  and  temperature 
profiles  do  not  peak  at  the  same  location  is  due  to  the 
growing  influence  of  the  viscous  forces  as  the  wall  is 
approached . 

The  above  argument  can  be  strengthened  somewhat  by 
appealing  again  to  the  equations.  Consider  the  differential 
equation  and  initial  conditions  for  the  dimensionless  velocity 
function  f(9). 


=  -  7M2  —  +  Re  —  -  (Ha)2  (<J>  -  f)  (105) 

dS2  3  g  g 


f(0)  =  1  ,  f'(0)  -  0  ,  g (0)  =  1 


Then  at  0  =  0 

f"(0)  =  -  7M2  +  Re  -  (Ha)2($  -  1)  (106) 

3 


For  the  moment  assume  that  the  temperature  function  g(0) 
either  remains  constant  or  decreases  with  9.  In  this  case 
if  the  parameters  Ha,  Re  and  M  are  so  chosen  that 


f"(0)  >  0 


i.e.  if 

-  7M2  +  Re  >  (Ha)2  (<t>  -  1)  (107) 

3 


then  the  function  f (0)  has  a  local  minimum  at  9  =  0.  Then 
as  9  increases  from  zero,  f  increases  and  all  terms  on  the 
right  side  of  Eq.  (105)  become  more  positive,  the  last  term 
becoming  less  negative.  Thus,  in  this  case,'  the  velocity 
function  f (9)  increases  monotonically  with  9,  growing  without 
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limit,  and  no  meaningful  solution  can  be  obtained.  Then 
if  f "(0)  >  0  meaningful  solutions  can  be  obtained  only 
if  the  temperature  function  g(9)  initially  increases  with  9 
and  increases  faster  than  f2.  Thus  in  order  to  have  a  bulge 
in  the  velocity  profile  there  must  also  be  a  bulge  in  the 
temperature  profile.  This  conclusion  supports  the  argument 
given  previously. 

Fig.  15  shows  the  velocity  gradient  at  the  wall  plotted 
against  $  for  the  three  values  of  Ha  used  in  this  study  and 
reflects  the  flattening  of  the  velocity  profiles  with  in¬ 
creasing  Ha  and  <J>.  This  flattening  of  the  velocity  profile 
and  the  corresponding  flattening  of  the  temperature  profile 
would  be  a  desirable  characteristic  for  wind  tunnel  accelerators. 
For  aerodynamic  testing  purposes  it  is  mandatory  that  the 
flow  approaching  the  model  be  quite  uniform. 

In  the  present  investigation  the  flattening  of  the 
profiles  is  due  not  only  to  the  increasing  ratio  of  magnetic 
forces  to  viscous  forces  as  represented  by  the  Hartmann 
number  but  also  to  the  corresponding  change  in  the  ratio 
of  inertia  forces  to  viscous  forces  as  represented  by  the 
Reynolds  number.  As  has  been  discussed  it  is  not  possible 
to  assess  the  influence  of  each  of  these  effects  separately. 

It  is  reasonable  to  suppose,  however,  that  with  larger 
values  of  the  Hartmann  number  than  used  in  this  study  both 
the  velocity  and  temperature  profiles  would  become  more 
uniform. 

Finally  a  comment  is  in  order  on  the  magnitudes  of  the 
values  of  Ha  and  Re.  In  the  formulation  of  this  problem 
these  parameters  were  defined  based  on  rQ,  the  radial 
distance  to  the  channel  entrance,  i.e. 


Ha  =  (-)2  Br 
0  ° 


Re  = 


p  U  r 
r  o  o  o 


A  more  conventional  procedure  would  be  to  define  these 
quantities  based  on  the  channel  height  at  rQ,  say  dQ,  i.e. 


«**  ■  <*>*  Vo 
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p  U  d 

Re*  = 


These  are  not  convenient  definitions  for  this  problem.  For 
small  values  of  the  divergence  angle  a 


d  fa  2ar 
o  o 


so  that 

Ha*  =  2aHa 
Re*  =  2aRe 


Then  for  a  =  1.5°  =  .02618  radians 


Ha*  =  .05236  Ha 
Re*  *=  .05236  Re 


Thus  the  values  of  the  Hartmann  number  and  Reynolds  number 
based  on  the  initial  channel  height  are  much  smaller  than 
the  values  based  on  the  radius  to  the  channel  entrance. 

For  small  values  of  Re*,  one  must  check  the  validity  of  a 
continuum  analysis  in  terms  of  the  Knudsen  number 


M 

Re* 


In  particular  for  Ha  —  100  and  $  —  1.2,  the  maximum  value 
of  M  and  minimum  value  of  Re*  used  in  this  investigation 
give  a  Knudsen  number  of  0.11.  This  value  suggests  that 
under  these  conditions,  the  continuum  analysis  is  perhaps 
being  pushed  to  the  limit  of  validity. 
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Vi 

VARIABLE  PROPERTIES 

6.1  FORMULATION  OF  THE  EQUATIONS 

It  was  shown  in  Section  III  that,  as  a  consequence  of 
the  separation  of  variables  technique  employed  in  this  study, 
the  electrical  conductivity  a,  the  viscosity  p.,  and  the 
thermal  conductivity  /c  cannot  vary  down  the  channel.  Equa¬ 
tions  have  been  formulated  and  solved  assuming  that  these 
transport  properties  are  constant  throughout  the  flow  field. 
Such  an  assumption  is  not  realistic  since  the  values  of 
these  properties  are  strongly  dependent  on  the  local  tempera¬ 
ture.  The  results  of  the  investigation  for  constant  properties 
indicate,  however,  that  the  temperature  variation  across  the 
channel  is  much  more  pronounced  than  the  temperature  variation 
down  the  channel,  i.e.  | ST/S0 | » | &T/Sr | .  It  was  felt,  then, 
that  useful  information  about  the  effect  of  the  variation 
of  the  transport  properties  on  the  velocity  and  temperature 
profiles  could  be  obtained  by  allowing  these  properties  to 
vary  as  a  function  of  temperature  across  the  channel  while 
ignoring  their  variation  down  the  channel. 

It  should  be  pointed  out  that  although  generally 
|  (3T/c30  |  »  |  t)T/c)r  j  ,  the  total  temperature  change  between  the 
entrance  and  the  exit  of  the  channel  may  be  quite  large. 

For  example,  if  L  is  the  channel  length  and  dQ  is  the  channel 
height  at  the  entrance,  then  the  ratio  of  centerline  tempera¬ 
ture  at  exit  to  the  centerline  temperature  at  entrance  is 


T  r  2  r  +  d  L/d^  „ 

-2.  =  (_2)  -  (_2 _ 2 _ 2) 

T  r  r 


(108) 


T 

e 


T 

o 


(1  +  2a  — )2 


(109) 


Then  for  a  =  1°  =  .01745  radians  and  L/d  =  10, 

o  ' 
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T 

—  =  1.82 


a  temperature  change  that  could  not  be  neglected  in 
evaluating  the  transport  properties.  Thus  the  equations 
that  will  be  developed  in  this  section  should  perhaps  be 
considered  valid  only  for  a  short  section  of  the  channel 
Ar  <  dQ. 

The  governing  differential  equations  will  first  be 
formulated  taking  into  account  only  the  variation  of  the 
electrical  conductivity.  Assume  a  a  -  T  relation  of  the 
form 


OJ 


(110) 


where  oc  and  Tc  are  reference  values  of  the  conductivity 
and  temperature,  say  values  along  the  channel  centerline. 
oc  is  assumed  to  be  constant.  Introducing  Eq.  (110)  into 
Eqs,  (34)  and  (36)  and  non-dimensionalizing  yields 


i  7m2  —  +  Re  —  -  (Ha)2g“(<J>  -  f)  (111) 

aep  3  g  g 


*L-&  -f  4g  =  2PrRef  -  (7  -  l)M2Pr  (Ha)2ga>(0  ~  f)2 
d62 


-  (7  -  l)M2Pr  (— )2  +  -  (37  +  1  )M2Prf  (112) 

de  3 


where  now 


Ha  =  (i)*  B0r0 
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The  initial  conditions  for  f  and  g  are 


f(0)  =  1,  f'(0)  -  0 

gCO)  =  1,  g'(0)  =  0  (113) 


Eqs.  (Ill)  and  (112)  together  with  initial  conditions  of 
Eq.  (113)  are  amenable  to  solution  by  the  same  technique 
as  used  for  the  constant  conductivity  model.  Figs.  44 
through  49  show  the  divergence  angle  a  and  the  ratio  of 
wall  temperature  to  centerline  temperature  g(a)  plotted  as 
functions  of  the  Reynolds  number  and  Mach  number  for 
o>  =  3/2  and  (Ha,4>)  =  (100, 1.6) , (150, 1.2) , (150, 1.6) . 

Equations  will  next  be  formulated  including  the  cross 
channel  variations  of  viscosity  and  thermal  conductivity 
according  to  a  square  root  law,  i.e. 


JL  =  JL 

U,  K 

^c  c 


(114) 


One  additional  assumption  is  made  that 


&  =  0 


(115) 


This  assumption  is  convenient  since  the  9  momentum  equation 
can  no  longer  be  directly  integrated  and  is  justified  by 
the  results  for  the  constant  property  model.  Then  intro¬ 
ducing  Eqs.  (110),  (114),  and  (115)  into  Eqs.  (34)  and 
(36)  and  non-dimensionalizing  gives 


d!f  =  Re  -11-  -  (M)(d£)  -  (Ha)2g^<*  -  f)  (116) 

d02  g3/S  2g  dg  d0 

+  4g  »  2RePr  ~4- —  (^&)2 

d@2  g^  2g  d0 
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-  (7  -  l)M2Pr (Ha)?gU)  *($  -  f)2 


-  (7  -  l)M2Pr 


+  £  z* 

dO  3 


(117) 


where 


Ha 


C  A 
(— )2 
M-r, 


B  r 
o  o 


Re 


PqVo 

P’c 


C  11 

Pr  =  -Elis. 


The  initial  conditions  on  f  and  g  are 


f(0)  -  1, 

f'(0)  «  0 

rH 

11 

O 

w 

to 

g'<0)  =  0 

(118) 


Again  Eqs.  (116)  and  (117)  with  initial  conditions  of 
Eq.  (118)  can  be  treated  by  the  technique  described  in 
Section  IV.  Figs.  50  through  55  show  a  and  g(a)  plotted 
as  functions  of  Re  and  M. 

6.2  DISCUSSION  OF  RESULTS 

The  results  of  this  portion  of  the  investigation  are 
presented  in  Figs.  20  through  25.  Figs.  20  through  22  show 
velocity  and  temperature  profiles  for  c  ~  T-5/2,  and  Figs. 

23  through  25  show  velocity  and  temperature  profiles  for 
o'  ~  Ts,/2,  ft  ~  (c  ^  T1/2 .  In  both  cases  the  divergence  angle 
is  approximate ly  1.25°  and  the  ratio  of  wrall  temperature  to 
centerline  temperature  is  approximately  0.2.  Figs.  17 
through  19  show  the  corresponding  velocity  and  temperature 
profiles  for  constant  properties.  Table  III  gives  the 
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comparative  values  of  the  temperature  and  velocity  gradients 
at  the  wall. 

The  most  conspicuous  feature  of  the  results  for  the 
variable  conductivity  model  is  the  accentuation  of  the 
bulges  in  both  the  velocity  and  temperature  profiles.  A 
comparison  of  Figs.  19  and  22  best  illustrates  this 
behavior. 

The  accentuation  of  the  temperature  bulge  can  be 
attributed  to  the  increased  Joule  heating  in  the  high 
temperature  region.  Using  Ohm's  law,  the  expression  for 
the  Joule  heating  can  be  written  as 


—  =  0Uo2Bo2($  -  f )2 


(119) 


Since  a  ~  T3/2  it  is  clear  that  this  term  will  be  large  in 
regions  of  high  temperatures. 

With  regard  to  the  velocity  bulge,  in  addition  to  the 
density  change  across  the  channel,  which  is  more  pronounced 
in  this  case,  there  is  also  an  increase  in  the  J  x  B  force 
in  the  high  temperature  region  as  a  result  of  the  conduc¬ 
tivity  variation. 

It  is  expected  that  for  larger  values  of  Ha  and  <t>, 
where  it  would  be  necessary  to  approach  the  limiting  values 
of  Re  and  M  in  order  to  maintain  a  constant,  both  the  velocity 
and  temperature  profiles  would  become  more  uniform. 

The  main  effects  of  the  variation  of  viscosity  and 
thermal  conductivity  are  to  depress  the  bulges  in  the 
velocity  and  temperature  profiles  and  to  steepen  the 
velocity  and  temperature  gradients  at  the  wall.  This  can 
be  seen  by  comparing  Figs.  25  and  26. 

In  the  region  of  high  temperature  the  corresponding 
high  value  of  thermal  conductivity  tends  to  promote  the 
diffusion  of  heat  to  the  adjacent  low  temperature  regions. 

Near  the  cool  wall  the  lower  thermal  conductivity  opposes 
the  transfer  of  heat. 

The  influence  of  the  viscosity  variation  is  analagous. 
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VII 

SUMMARY 


A  solution  has  been  obtained  to  the  magnetogas dynamic 
equations  of  motion  including  the  effects  of  compressibility, 
viscosity,  thermal  conductivity  and  crossed  electric  and 
magnetic  fields.  The  solution  obtained  is  an  exact  solution 
in  the  sense  that  no  terms  have  been  dropped  from  the 
governing  equations  and  boundary  layer  approximations  have 
not  been  employed. 

The  particular  flow  model  considered  was  the  two- 
dimensional,  fully  developed  flow  between  diverging  plane 
walls.  Self-similar  solutions  to  the  governing  equations 
were  sought.  It  was  determined  that  such  solutions  could 
be  obtained  under  certain  conditions;  the  resulting  flow 
being  an  accelerating  flow,  the  velocity  increasing  linearly 
with  distance  down  the  channel,  in  which  the  pressure  and 
Mach  number  are  constant  along  streamlines.  For  this  special 
type  of  flow,  with  the  assumption  of  constant  transport 
properties,  the  partial  differential  equations  of  motion 
reduce  to  a  pair  of  coupled,  non-linear,  ordinary  differ¬ 
ential  equations,  Eqs.  (45)  and  (46)  with  boundary  conditions 
given  by  Eq.  (47) . 

It  was  shown  that  the  boundary  value  problem  defined 
by  Eqs.  (45)  through  (47)  is  generally  overspecified.  That 
is,  if  the  divergence  angle  a  and  the  ratio  of  wall  temper¬ 
ature  g(a)  are  considered  fixed,  and  two  of  the  four  para¬ 
meters  Ha,  Re  and  M  are  specified,  then  the  remaining  two 
parameters  are  uniquely  determined.  Thus,  in  a  broad  sense, 
Eqs,  (45)  through  (47)  constitute  an  eigenvalue  problem. 

If  the  required  relationship  among  the  parameters  is  not 
satisfied  then  a  self-similar  solution  to  the  governing 
equations  does  not  exist.  The  technique  used  for  the  solution 
of  this  problem  has  been  amply  discussed  in  Section  IV. 

The  results  of  the  investigation  were  presented  in  the 
form  of  velocity  and  temperature  profiles  for  several  values 
of  the  electromagnetic  parameters  Ha  and  <£>.  The  values  of 
the  gas  dynamic  parameters,  the  eigenvalues,  were  selected 
so  as  to  maintain  the  divergence  angle  and  the  ratio  of  wall 
temperature  to  centerline  temperature  constant.  The  most 
conspicuous  feature  of  the  results  is  the  prominent  bulge 
that  occurs  in  the  temperature  distribution  for  certain 
values  of  the  parameters.  This  behavior  becomes  less  pro¬ 
nounced  for  large  values  of  Ha  and  4>  where,  incidentally, 
practical  devices  would  be  likely  to  operate. 
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Although  it  was  not  possible,  in  this  study,  to  include 
the  variation  of  the  transport  properties  a,  jj,  and  k  along 
the  channel,  it  was  shown  that  the  effect  of  the  variation 
of  these  properties  across  the  channel  is  to  accentuate  the 
bulges  in  both  the  velocity  and  temperature  profiles. 

Finally,  a  few  comments  are  in  order  on  the  significance 
of  the  present  study  with  regard  to  wind  tunnel  accelerators. 
Certainly  the  mathematical  model  used  in  this  investigation 
is  too  crude  to  be  utilized  for  detailed  design  of  a  prac¬ 
tical  accelerator.  Nonetheless  there  are  certain  qualitative 
features  of  the  results,  such  as  the  bulge  in  the  temperature 
profiles,  which  are  likely  to  be  encountered  in  a  practical 
device.  Also  there  are  several  characteristics  of  the  flow 
situation  considered  in  this  study  which  might  be  desirable 
for  wind  tunnel  accelerators.  For  example,  although  a  large 
portion  of  the  energy  added  to  the  plasma  is  dissipated  in 
Joule  heating  thereby,  increasing  the  enthalpy  of  the  gas, 
this  enthalpy  would  be  recoverable  as  directed  kinetic  energy 
in  the  post  accelerator  nozzle.  In  particular  this  expansion 
should  pose  no  difficulty  since  no  pressure  drop  occurs  in 
the  accelerator  itself.  Furthermore,  despite  the  fact  that 
the  flow  is  considered  fully  viscous,  for  large  values  of 
Ha  and  0  a  rather  well  defined  core  region  does  exist  where 
the  velocity  and  temperature  profiles  are  relatively  flat. 

From  these  considerations,  it  would  appear  desirable 
to  extend  the  present  investigation,  using  a  more  realistic 
model.  One  might,  for  example,  relax  the  constraint  of 
separation  of  variables  technique  employed  here.  This 
technique,  despite  its  undeniable  virtue  of  simplifying  the 
problem,  imposes  rather  severe  restrictions  on  the  solution 
obtained.  Then,  rather  than  assuming  a  particular  form  for 
the  solution  as  was  done  in  Section  III,  one  would  specify 
initial  conditions  at  the  channel  entrance  rQ.  This  pro¬ 
cedure  would  require  solution  of  a  pair  of  coupled,  non¬ 
linear,  partial  differential  equations  with  boundary  con¬ 
ditions,  a  very  difficult  but  not  impossible  task.  It  is 
likely  that  some  of  the  well  known  techniques  developed 
in  boundary  layer  theory  could  be  extended  to  this  type  of 
channel  flow  problem.  Such  a  solution  would  give  detailed 
information  for  a  realistic  MHD  channel  flow  to  a  degree  not 
presently  available  in  the  open  literature. 
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Post -Accelerator 
Nozzle 


Fig.  1  MHD  Powered  Wind  Tunnel 
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Ha  =  100 

a  =  1.503 

<*>  =  1.2 

g(a)  =  .200 

M  =  1.87 

Re  -  350 

Fig.  6  Velocity  and  Temperature  Profiles  for  Ha  =  100,  =  1.2 
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Ha  =100  a  =  1.503° 

<t  =  1.6  g  (a)  =  .209 

M  =  2.24 
Re  =  5150 


Fig.  7  Velocity  and  Temperature  Profiles  for  Ho  =  100,  0  =  1.6 
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Ha  =  100  a  =  1.502° 

0  =  2.0  g (a)  =  .202 

M  =  2.02 

Re  =  9700 


Fig.  8  Velocity  and  Temperature  Profiles  for  Ho  *=  100,  $  =  2.0 
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Ha  =150  a  -  1.500° 

<t>  -  1.2  g(a)  =  .207 

M  =  3.0 
Re  =  4100 
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Ha  =150  a  =  1.502° 

<t>  =  1.6  g(a)  =  .202 

M  =  2.704 
Re  =  14010 


Fig.  10  Velocity  ond  Temperature  Profiles  for  Ha  =  150,  =»  1.6 
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Ha  =  150  a  =  1.500° 

<t>  -  2.0  g(a)  =  .202 

M  =  2.209 

Re  =  23045 


0  0.2  0.4  0.6  0.8  1.0 


Fig.  II  Velocity  and  Temperature  Profile*  for  Ho  =  150,  4>  =  2.0 
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Ha  =  200  a  =  1.497° 

<t>  =  1.2  g(a)  =  .  191 

M  =  3.75 

Re  =  8450 

f (n>  g(n) 


0  0.2  0.4  0.6  0.8  1.0 


n 

Fig.  12  Velocity  and  Temperature  Profiles  for  Ha  =  200,  <t>  =  1.2 
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Ha  =  200  a  =  1.501° 

=  1.6  g  (a)  =  .  190 

M  =  2.86 

Re  =  24600 


0  0.2  0.4  0.6  0.8  1.0 


n 

Fig.  13  Velocity  ond  Temperature  Profiles  for  Ho  =  200,  <t>  =  1-6 
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Ha  -  200 

a  =  1.504 

O  =  2.0 

g (a)  =  .202 

M  =  2.2376 

Re  =  40250 

f(il)  g(t}) 


Fig.  14  Velocity  and  Temperature  Profiles  for  Ho  =  200,  *=  2.0 
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Ha  =  100 
<D  =  1.6 
M  =  1.95 
Re  =  3800 


a  =  1.252° 
g(a)  =  .206 


Fig.  17  Velocity  and  Temperature  Profiles  for  Ho  =  100,  $  “  1.6 
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Ha  =150  a  =  1.250° 

<t>  =  1.2  g(a)  =  .197 

M  =  2.515 
Re  =  3075 


0  0.2  0.4  0.6  0.8  1.0 


TJ 

Fig.  18  Velocity  and  Temperature  Profiles  for  Ho  =  150,  $  ”  1.2 
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Ha  =  150 
4>  -  1.6 
M  =  2.525 
Re  =  13490 


a  =  1.249 
g<a)  =  .199 
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Ha  =  100  a  =  1.249° 

<J>  =*  1.6  g  (a)  =  .200 

M  =  1.96 
Re  =  3925 


71 


2 


AEDC-T  R-67-80 


Ha  =150  a  =  1.245° 

<t>  =  1.6  g(a)  =  .199 

M  =  2.815 
Re  ^  17500 
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Ha  =100  a  =  1.249° 

<t>  =  1.6  g(a)  =  2.02 

U  =  2.05 
Re  =  4550 
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Ha  =150  a  =  1.251° 

<b  =  1.2  g  (ct)  =  .197 

M  =  2.71 
Re  =  4200 


Fig.  24  Velocity  end  Temperature  Profiles  for  Ha  =  150,  <t>  =  1.2, 

X5'2  x l-"2 

a  ~  T  ,  jJ.  ~  k  ~  i 
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Fig.  27  Variation  of  g(a)  with  Re  and  M  for  Ha  =  100,  O  =  1.2 
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Fig.  35  Variation  of  g(a)  with  Re  and  M  for  Ha 
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Fig.  42  Variation  of  a  with  Re  and  M  for  Ha  =  200,  4*  =  2.0 
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Fig.  52  Voriationof  a  with  Re  and  M  lor  Ha 
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TABLE  I 

LIMITING  VALUES  OF  Qi 


M  AND  Re 


Ha 

4> 

a 

min 

Mlimit 

Ren  .  . 

limit 

100 

1.2 

1.420° 

4.9894 

1906.6 

100 

1.6 

.938° 

2.8857 

5968.9 

100 

2.0 

.756° 

2.2357 

9981.4 

150 

1.2 

.946° 

4.9961 

4406.6 

150 

1.6 

.625° 

2.8863 

8068.9 

150 

2.0 

.504° 

2.2359 

22481.4 

200 

1.2 

.710° 

4.9973 

7906.6 

200 

1.6 

.469° 

2.8865 

23968.9 

200 

2.0 

.378° 

2.2360 

39981.4 

TABLE  II 

SUMMARY  OF  RESULTS  FOR  CONSTANT  PROPERTIES 


Ha 

Re 

M 

f '( a) 

g'(a) 

P*(0) 

100 

1.2 

350 

1.87 

-  117 

-  186 

.0372 

100 

1.6 

5150 

2.24 

-  138 

-  347 

.00361 

100 

2.0 

9700 

2.Q2 

-  156 

-  386 

.00157 

150 

1.2 

4100 

3.00 

-  172 

-  629 

.00819 

150 

1.6 

14010 

2.704 

-  209 

-  776 

.00194 

150 

2.0 

23045 

2.209 

-  238 

-  711 

.00079 

200 

1.2 

8450 

3.750 

-  231 

-  1333 

.00621 

200 

1,6 

24600 

2.860 

-  281 

-  1180 

.00124 

200 

2.0 

40250 

2.2376 

-  318 

-  980 

.00046 

TABLE  111 

COMPARISON  OF  i'[a)  AND  g'{«)  FOR  CONSTANT  AND  VARIABLE  PROPERTIES 


Constant 


Properties 

Ha 

4> 
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1.6 

-  207 

-  666 
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T3/2 

o 

p  K 

r\j  T3/2 

<v  T3/2 

(a) 

s'  (a) 

f  '  (a) 

g'(ct) 

110 

-  204 

-  251 

-  499 

181 

-  576 

-  392 

-  1141 

322 

-  1322 

-  662 

-  2434 
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